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ABSTRACT 


An  experimental  investigation  was  undertaken  in  the 
low  density  wind  tunnel  to  test  the  validity  of  the  boundary 
layer  equations  when  used  with  "slip"  and  "temperature  jump" 
boundary  conditions  in  the  transition  regime. 

The  boundary  layer  equations  were  solved  using  a 
numerical  technique  developed  by  D.  Clutter  and  A.M.O.  Smith, 
using  as  boundary  conditions  the  classical  slip  and  tempera¬ 
ture  jump  boundary  conditions,  and  also  by  using  experi¬ 
mentally  determined  slip  and  temperature  jump  results.  The 
slip  theory  predicts  lower  values  of  slip  than  were  experi¬ 
mentally  obtained,  and  the  boundary  layer  profiles  begin  to 
correlate  with  experiment  after  approximately  one  inch  from 
the  leading  edge,  where  the  local  Reynolds  number  is  of  the 
order  of  11. 

Inclusion  of  temperature  jump  in  the  theory  indicates 
that  its  effect  on  the  slip  results  is  small  at  least  for 
Knudsen  number  less  than  0.1. 

Changing  the  initial  value  of  slip  velocity  at  the 
leading  edge  does  not  change  the  slip  profile  appreciably 
after  three  or  four  mean  free  paths  from  this  point. 
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P,P 
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X 

a 
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G 
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m  mass  of  gas  molecule 

’  denotes  differentiation  with  respect  to  0 
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Other  symbols  are  defined  where  they  appear 


CHAPTER  I 


INTRODUCTION 


An  analytical  solution  of  the  equations  governing  a 
physical  system  is  facilitated  by  a  mathematical  model  which 
describes  the  behavior  of  such  a  system.  The  relationships 
describing  such  a  mathematical  model  are  known  as  constitu¬ 
tive  equations.  The  constitutive  equations,  which  are  neces¬ 
sary  for  a  solution  to  the  problem,  incorporate  the  material 
response  in  the  formulation  of  the  problem.  Eringen  [1]* 
notes  that  "The  character  of  the  material  is  brought  into  the 
formulation  through  appropriate  constitutive  equations  for 
each  material  with  the  constitutive  variables  being  restricted 
in  their  regions  of  definitions". 

Certain  material  properties,  such  as  elasticity,  visco¬ 
sity,  and  plasticity,  appear  as  parameters  in  the  constitutive 
equations;  these  properties  are  termed  by  Reiner  [2]  as  funda¬ 
mental;  and  other  properties,  such  as  visco-elasticity  and 
visco-plasticity,  as  complex  properties.  The  fundamental 
material  properties  may  be  associated  with  the  manner  in 
which  mechanical  energy  is  dissipated  by  the  material.  In 
the  case  of  gases  at  low  pressure  the  linear  Newtonian  model 
for  the  relationship  between  shear  stress  and  shear  strain 


*  [  ]  denotes  reference. 
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rate  is  appropriate. 

The  constitutive  equation  for  simple  shear  of  a  New¬ 
tonian  fluid  is , 


Txy  =  2\i  dxy 

where  Txy  is  the  shear  stress  and  dxy  is  the  shearing 
strain  rate.  The  viscosity  coefficient,  y  ,  depends  only 
on  temperature  and  pressure,  and  not  on  the  rate  of  deforma¬ 
tion.  The  Newtonian  fluid  dissipates  energy  by  the  trans¬ 
fer  of  momentum  on  a  molecular  scale. 

A  theoretical  investigation  of  flow  problems  involves 
the  solution  of  the  equations  of  motion,  along  with  the 
energy  equation. 

The  equations  of  motion  are: 


9o  xx  ,  8  T  xy  ,  9  t  xz 

9  x  9y  9  z 


pax  » 


+  +  +  pfy  =  pay  , 


3txz  +  +  3ozz  +  pf 

3z  3y  3z 


Paz 


where  fx,  f  y  ,  fz,  and  ax,  ay ,  and  az  are  the  x,  y,  and 
z  components  of  the  body  force  and  acceleration  vector 
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respectively.  Furthermore,  if  the  material  under  considera¬ 
tion  is  a  Newtonian  fluid,  then  the  governing  equations  may 
be  further  reduced  to  the  Navier-S tokes  equations,  which  in 
vector  notation  are  given  by  [3] 

P  Dt  =  p7  “  7P  _  f  V^V-V)  +  2  (Vy)  .  VV 

+  (Vy)  x  (VxV)  +  y[V2V  +  V(V.V)] 

where  V  is  the  velocity  vector,  f  is  the  body  force  per 
unit  mass,  and  P  is  the  static  pressure.  The  Navier- 
Stokes  equations  are  not  readily  solvable  in  their  general 
form.  However,  in  some  important  special  cases  these  equa¬ 
tions  reduce  to  a  solvable  form. 

The  Knudsen  number,  Kn  ,  is  a  non-dimensional  para¬ 
meter  which  is  important  in  rarefied  gas  flows,  although  not 
derivable  from  the  Navier-S tokes  equations. 

The  Knudsen  number  which  is  defined  by  [4]  as, 


where  A  is  the  mean  free  path  and  D  a  characteristic 
length  which  will  be  taken  to  be  the  plate  length  in  this 
study,  distinguishes  between  flow  regimes  within  the  laminar 
region. 

Depending  on  the  relative  magnitude  of  the  Knudsen  number, 


< 
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the  laminar  flow  regime  may  be  divided  [5]  into  the  follow¬ 
ing  regimes: 

(a)  the  continuum  regime,  which  is  defined  when  Knudsen 
number  is  less  than  0.01, 

(b)  the  free  molecular  regime,  which  occurs  when  Kn  >  1, 
and 

(c)  the  transition  regime,  which  is  the  regime  in  which  a 
transition  occurs  between  the  continuum  and  free  mole¬ 
cular  regime. 

The  regime  of  interest  in  the  study  presented  in  this 
thesis  is  the  transition  regime;  experimental  investigations 
were  carried  out  at  Kn  *  0.11  and  Kn  *  0.06. 

Since  this  thesis  primarily  uses  the  Navier-S tokes 
equations  in  the  transition  regime,  it  is  important  that 
the  validity  of  the  use  of  these  equations  in  this  regime 
is  investigated.  The  most  important  point  to  be  considered 
is  the  boundary  conditions.  The  "no-slip"  boundary  condition 
is  usually  assumed  for  Newtonian  viscous  flow  and  its  range 
of  validity  is  defined  by  Knudsen  number  less  than  0.01. 

The  no-slip  boundary  condition,  as  assumed  in  classi¬ 
cal  aerodynamics,  is  clearly  inconsistent  with  the  observed 
behavior  [6]  of  the  flow  in  the  transition  flow  regime;  that 
is  for  Knudsen  number  in  the  range 


0.01  <  Kn  <  1  . 
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The  phenomenon  of  slip  was  first  observed  by  Kundt  and 
Warburg  [7].  These  investigators  found  that  the  damping  of 
a  vibrating  disk  by  a  surrounding  gas  decreased  at  low  pres¬ 
sures.  They  attributed  this  to  slip  of  the  gas  at  the 
boundaries.  It  was  not  until  four  years  later  that  a  theore¬ 
tical  expression  describing  this  phenomenon  was  obtained  by 
Maxwell  [8].  Maxwell  concluded  that  the  amount  of  slip  at 
the  wall  is  a  function  of  the  velocity  gradient  at  that  point, 
and  is  given  by 


Similar  conclusions  were  obtained,  but  in  a  simplier 
manner,  by  Millikan  [9]  and  reproduced  by  Kennard  [10].  Not¬ 
withstanding  some  slight  modifications  in  the  slip  coefficient, 
this  theory  is  still  widely  used  today.  It  has  been  shown 
experimentally  [11]  that  the  accommodation  coefficent  must 
depend  on  the  physical  condition  of  the  surface  and  on  the 
gas  under  consideration.  A  theoretical  expression  relating 
these  quantities  has  not  been  found  as  yet,  consequently  the 
slip  coefficients  used  in  any  theoretical  investigation  are 
based  on  experimental  results.  The  slip  coefficients  adopted 


in  this  investigation  are 
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a i  =  2.333  , 
a  x  =  1.857  , 

and 

ax  =  1.222  , 

which  correspond  to  the  following  tangential  momentum  accom¬ 
modation  coefficients 


and 


C  =0.6  , 

a  =  0.7  , 

a  =  0.9  respectively. 


These  values  cover  the  range  of  values  obtained  experiment¬ 
ally  by  Merlic  [12], 

So  far  attention  has  been  focused  on  the  velocity  boun¬ 
dary  condition.  The  other  boundary  condition  to  be  considered  is 
on  the  gas  temperature  at  the  surface.  In  the  continuum 
regime  the  gas  temperature  at  the  wall  is  equal  to  the  wall 
temperature.  However,  in  the  transition  regime  the  gas  tem¬ 
perature  at  the  surface  can  differ  from  the  wall  temperature, 
giving  a  "temperature  jump"  boundary  condition  [13].  The 
presence  of  such  a  temperature  jump  was  first  assumed  by 
Poisson,  and  later  confirmed  experimentally  by  Smoluchowski 
[14].  A  theory  describing  the  temperature  jump,  similar  to 
the  slip  theory  derived  by  Maxwell,  was  derived  by  Poisson. 


. 


*1-8 


7 


Poisson’s  temperature  jump  boundary  condition  is, 


» 


where  is  the  temperature  of  the  gas  that  would  exist  if 

the  temperature  gradient  along  the  outward-drawn  normal  to 
the  wall,  Al.  ,  continued  without  change  up  to  the  wall  itself 

9y 

[15],  and  g'  represents  the  "temperature  jump  distance". 

Again  it  is  clear  that  the  temperature  jump  is  a  func¬ 
tion  of  an  accommodation  coefficient,  which  is  referred  to  as 
the  thermal  accommodation  coefficient. 

It  is  worth  noting  that  the  thermal  accommodation  co¬ 
efficient  and  the  tangential  accommodation  coefficient  are 
not  necessarily  equal.  The  temperature  jump  coefficient  used 
in  this  thesis  is, 

a?  =  1.875  , 

which  is  based  on  the  thermal  accommodation  coefficient 

a  =  1 . 0  , 

and  is  the  coefficient  usually  assumed  for  diffuse  reflection. 

The  possibility  of  using  the  Navier-S tokes  equations 
subject  to  the  slip  boundary  condition  as  a  model  to  study 
viscous  flow  problems  in  the  transition  regime  was  investi- 


. 

'  f 
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gated  by  Patterson  [16]  and  by  Street  [17].  Street  con¬ 
cluded  that  the  slip  condition  is  compatible  with  the  Navier- 
Stokes  equations,  and  furthermore  that  it  may  be  possible 
to  study  viscous  flow  problems  in  the  transition  regime  us¬ 
ing  the  boundary  layer  equations,  which  take  the  first  order 
slip  effects  into  account. 

The  objective  of  this  work  is  to  experimentally  in¬ 
vestigate  the  range  of  validity  of  the  boundary  layer  equations 
and  transition  boundary  conditions. 


CHAPTER  II 


THEORETICAL  ANALYSIS 


2 . 1  Order  of  Magnitude  Analysis 

The  solution  of  many  flow  problems  can  be  facilitated 
by  an  order  of  magnitude  analysis  of  the  continuity,  energy, 
and  momentum  equations.  The  purpose  is  to  retain  those 
terms  which  are  essential  to  describe  the  flow. 

The  governing  equations  for  a  two  dimensional  steady 
flow  in  rectangular  Cartesian  coordinates  (see  reference 
[18]  or  [19]  for  derivation)  become: 

Cont inui ty : 


3(pU)  ,  3(PV)  _  o 
3X  3Y 


Longitudinal  Momentum 


3 P  ,  3 (-2/3  u)  . 
3X  3X 


' 


10 


Lateral  Momentum 


P(UH  +  vw> 


9P  ,  3 (-2/ 3  u)  ,3U  3V^ 

3Y  3Y  ax  dY 


+ 


(-  i 

v  3 


U) 


_a_ 

8Y 


(i£  +  av}  ajj  (av  au} 

V3X  3Y;  8x  ^3X  3Y; 


a2 


V 


+  2  +  y  ( 

3Y  3Y  ax2 


111)  +  ul  <M  +  - 

9y2;  m3Y  V3X  3Y; 


PG 


Energy : 


2.2-1 


p  C  p  ( U-j~  +  v*~)  +  u~-  +  v|£ 

H  H  ax  3Y  ax  3Y 


.  ra2T  a2Tx  ak  3t  ak  a t 

ax2  3y2  ax  ax  3y  3y 


+ 


(3U  aV)2  2u[(M)2 

^ax  ax;  yuax; 


+ 


(av)2]  (au 

V3Y;  J  MV3Y 


+ 


avx  2 
ax' 


The  boundary  conditions  at  the  wall  will  be  taken  to 
be  the  Maxwell  slip  condition,  and  for  convenience,  the 
temperature  of  the  gas  at  the  wall  will  be  assumed  constant. 


Therefore,  at  Y  *  0 


11 


U(x,0) 


X  a, 


au 

9Y 


where  a-^  is  a  constant. 


T (X, 0)  -  Tw  2.1-2 

V(X,0)  ■  0  (no  blowing  or  suction) 

P (X, o)  =  Pw 


At  Y  =  00  : 


U 

T 

V 

P 


=  U 

=  T 

-  0 

-  P 


00 


00 


The  terms  in  the  equations  are  nondimensionalized  by 
dividing  all  the  variables  by  characteristic  quantities, 
denoted  by  subscript  C. 


Pc_Uc 

X 


<*£  *  ■£>  - 


0 


2 

£] 


9u  ,  rPcvc^Ci  9u 

“U3l  +  yT  “V3? 


9P 

9x 


+ 


,4  9M"  9u  ,  4  9 2 u v 

(3  TZ  +  3  M 


' 


.  - .  • 
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+  r  ^ c V c  i  ,  2  9M"  9_v  ,  M"  92v  9M'  3V 

XCYC  3  3x  9y  3  9x9y  9y  9x 


+ 


y.u 


cuc 


,9M" 

°y 


9  u 

9y 


+  M 


92Ux 

9y2 


Pcucvc 

Xc 


c 


c 


2 

-] 


9_P 

9y 


,  rycUCi  /  I  3*T  9u  M"  92 u  9M"  9u>. 

lYcXcJ  v  3  9y  9x  3  9y9x  9x  9y' 


,  r^c^Ci  /4  9M^  9v  ,  4  w^92vn 

+  L - tJ  \rr  7T—  -r—  +  -r  M  - =0 

Y  2  3  9y  9y  3  9y2 


2.1-3 


+  [ilsZ|]  (|T  |2  +  M'l^)  -  [PcO]» 


X, 


2  9x  9x  9x2 


PcCpUc0Cl  90  .  cSVC®C!  90  rUC^Ci  9P  rVCPC1 

[ - £ - J  WUtt  +  [ - £ - J  (a)  V-r— *■  -  [— ]  u-r —  -  [— ]  V' 


9x 


9y  Xc  9x 


9P 

9y 


[lc£c,  (Ki±!  +  |JL  |4)  +  [l£Zf]  (Kll|  +  M  |1) 

xc2  9x2  9x  9x  Y c 2  9y 2  9y  9y 


KJ, 


+  [Mfi]  iS'  (|^)2  +  [^|i]  is'  (|^)2 

x  2  3  9x  y  2  3  9y 
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+ 


ycvcuc 

XCYC 


]  (- 


4M^  3u  3v  ,  0Vf^3u  3vn 

~3  37  77  M  3?  3^ 


+ 


_c 

2 


2 


]  M" 


(lii)2  + 
3y 


2 

c  Y  c  .] 


y~v 


M' 


(3v)2 

dx 


[Uc]  u (x , 0) 


where 


K  =  K/Kc 

m"  =  y /yc 

a)  =  p/pc 

e  =  t  -  T  oo 
4>  =  e/ec 


As  s  ume 


Uc  =  Uoo 

Xc  =  L  (the  plate  length) 

Pc  -  poo 

8C  *  Tw  -  Too  =  ew 


. 
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Uc  *  Uoo 


If  it  is  assumed  that  in  the  laminar  boundary  layer 
the  order  of  magnitude  of  the  viscous  forces  is  equal  to  the 
order  of  magnitude  of  the  inertia  forces,  this  implies. 


or  Y 


j  PooUqO 


- t-tt*  where  ReT  is  the  Reynolds  number 

(Re,)1/2  L 


based  on  plate  length. 

However,  from  the  slip  boundary  condition, 


1 


Therefore  let, 


Y 


c 


KnL  + 


(ReL)1/2 


0  <  Kn  <  1 


All  terms  in  the  continuity  equation  are  taken  to  be 


of  equal  order.  Then 


.riugnsl  ai*Iq  If*  fc*l 
,,  ari3  .iirotl  ,*•*•*•« 
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or 


P  U 
K00  00 


PcoVc 

=  o  (— — — — ) 


v„  = 


UooY, 


=  Uoo  (Kn  + 


(ReL> 


1/2 


Substituting  the  values  for  the  characteristic  quanti¬ 
ties  into  the  equations  and  re-arranging 


Hr  +  «¥■  + 

9x  9x 


U> 


9  v 

9y 


+  u|£  = 

3y 


o 


9  u  .  9  u 

“u3l  +  “V3? 


-  [Eu]  f£  +  [-^]  <±  |S' 

9x  Re^  3  9x  9x  3  ^2 


+  |-  1  ]  2  9M"  9v  +  M"  9 2v  +  9M"  9v^ 


Re, 


3  9x  9y  3  9x9y  9y  9x' 


Re- 


+  i 


(Kn  +  (l-Kn)/(ReT)1/2)2J  ^  ^ 


,  ,9M^  9u  ^92us 

]  +  * - T) 


9y 


0)U 


9  v 
9x 


+  03V 


9  v 

9y 


[ 


Eu 


(Kn  +  (l-Kn)/(ReL) 


] 


1/  2  }  2 


9P 

3y 
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+ 


1 


[ 


Re- 


(Kn  +  ( 1-Kn )  /  ( ReT  )1/2)2 


2  8M*  3 u 

3  3y  3 x 


+ 


M’  32  u  , 

3  3  y  3  x 


3M’  3u>. 
3x  3y 


1 

,  r _ ReL _ i  ,_4  3M*  3v  ,  4_M' _3_^_v, 

(Kn  +  (1-Kn) / (ReL) 1/2) 2  3  9y  9y  3  3y2 


Re- 


/  3_M*  _3_v 
3x  3x 


H^f) 

3  x2 


[-]oj 

lF 


2 . 1-4 


03  U' 


3  (f) 
3  x 


+  0)V 


ii 

3y 


-  [E 


EJ 


,  i p 

(u3x 


3P. 
v — ) 


3y 


=  [ 


P  r  Re- 


]  (K 


3  2  4> 

3x2 


3K  3(f) 
+  3x  3x' 


+  [ 


(Pr x ReL ) 


- -  r— - -]  (K-^  +  |^) 

(Kn  +  (1-Kn) / (ReT ) 1/2) 2  3y2 


,  r  E  w  AM'  iH.  AX  +  lX) 
+  [R^]  ("  T  37  37  +  2M  3y  3x; 


+  [■ 


Re. 


(Kn  +  (l-Kn)/(ReL)  ) 


_ 1  m'  (— ) 

1/2, 2J  M  V3y; 


*6 
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where 


+  [Kn  +  (l-Kn)/(ReL)1/2]2  (|i) 2 


Eckert  number,  E  = 


U 


00 


Cp  6W 


Euler  number,  Eu  = 


00 


P  00  U 


00 


u 


Froude  number,  F  = 


00 


LG 


c  y. 


Prandtl  number,  Pr  =  —2. 


K 


oo 


Consider  the  flow  where  E  £  0(1),  F  >  0(1),  Pr  =  0(1), 
and  0.5  <_  M  _<  1 ;  the  bounds  on  Mach  number  implying  that 
Eu  =  0(1) . 

From  the  " y " -momen t um  equation 


9P 

9y 


[Kn  +  (l-Kn)/(ReL)1/2]2 

0( - - - - - )  <  0(1). 


and 


ajp 

ay 


Since  P  is  assumed  to  be  constant  in  the  free  stream 
r)  P 

~  0 ,  then  tt —  ~  0  . 


£ 


U 


' 

00  s  .  ?<i 

ISC  1  u 

luhnsi' 

■ 

. 
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For  the  conditions  stated  above  the  dimensional  equa¬ 
tions  governing  the  flow  become: 


!x(pu)  +  l7(pv)  =  0 


2.1-5 


p[uM  +  vH]  =  1_  [yM] 

Pl  9X  3YJ  9Y  lP8YJ 


2.1-6 


p[U 


3h 

8X 


V— ] 
3YJ 


y( 


9  IK 
9y' 


(K 


2.1-7 


where  Cp 


dh 

dT 


The  boundary  conditions  at  the  wall  are: 


U(X,0)  =  Uw  (slip) 

T(X,0)  =  T ^  (temperature  jump) 

V(X,0)  =  0 

outside  the  boundary  layer, 

U  =  Uoo 
T  =  Too 
V  =  0 

and  at  X  =  0 

U  =  UM 
T  =  Too 


2.1-8 


■ 
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The  solution  to  the  problem  has  been  attempted  by 
Maslen  [20]  who  perturbated  around  a  similar  solution  in 
order  to  take  into  account  slip  and  temperature  jump  boun¬ 
dary  conditions.  This  approach  can  be  improved  upon  by 
numerically  solving  the  above  system. 


2 . 2  Numerical  Analysis 

Since  a  closed  form  solution  to  the  system  in  Section 
2.1  did  not  seem  possible  a  numerical  solution  of  the  equations 
was  attempted. 

Introducing  total  enthalpy, 


H 


» 


the  energy  equation  becomes, 


(U 


9H 

9X 


VM)  = 

9y' 


_9_ 

9Y 


(-U- 

Pr 


9H 

9Y 


+ 


uH) 

9y' 


2.2-1 


Following  the  numerical  method  presented  by  D.  Clutter  and 
A.M.O.  Smith  [21],  the  equations  are  transformed  using  a  modi¬ 
fied  Howarth-Dorodnitsyn  transformation  to  a  more  convenient 
coordinate  system.  The  transformation  has  the  advantage  that 
for  this  problem  the  singularity  at  x  =  0  is  removed  and  the 


-±bo*  a  aolao  baoio^araTd  9ia  enoi3aupa  adl  ,1  ]  d*  •  • 

■ 
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boundary  layer  thickness  in  the  transformed  plane  does  not 
grow  as  rapidly  along  the  plate  as  in  the  x,  y  coordinate 
sy s  tem . 

The  Howarth-Dorodni tsyn  transformation  is 


nH  =  /y  p  dy 

0 

£  =  x 

Clutter  and  Smith  have  introduced  a  modified  trans¬ 
formation  which  reduces  to  a  form  studied  in  previous  work 
on  incompressible  flow, 

Thus  , 


n 


u 


( 


Poo  y  oo  X 


»_)l/2 


P  dy 


X  =  X 


A  stream  function  is  defined  such  that 


pu 


3y 


pv 


8^ 

8x 


If  we  define  another  stream  function  f  such  that, 


or  I  mQ 
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if  _  f.  .  _u_ 

3n  u*, 


the  relationship  between  and  f  becomes 


*  =  -'PooUcoX  f  • 


Introducing  the  non-dimensional  variables,  g  =  and 

Hr 


l00 


c  = 


py 


Poo  Poo 


,  the  transformed  equations  become; 


3n 


[cf" 


+ 


ff " 

2 


x[f  1 


1 1 


M]  = 

8  x J 


0 


2 . 2-2 


3n 


[£  8'  +  fe  C(1-P7> 


f'f" 


la' 

2 


2 . 2-3 


+  x[f' 


All  primes  denote  differentiation  with  respect  to  r| . 

Round  off  errors  in  the  computer  program  can  be  reduced 
[21]  by  introducing  the  following  variables: 


,p  oJ  39qei)i  •  t  k\j  1C  iSsilne-isl  j  9-jc:-s  n  *q 
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<J> 

=  f-n 

i — i 

i 

00 

II 

CD 

=  f’-i 

0’  =  g' 

2.2-4 

4>” 

=  f" 

0”  =  gM 

<j>  "  ’  =  f  »» ' 

With  these  new  variables,  equations  2.2-2  and  2.2-3 
b  e  come : 


_3_ 

3n 


[Ccj)M 


( 4>+n )  <$> " 
2 


2 . 2-5 


+  x[(f+i)  M.' 

3x 


Ml 

3  x  J 


_3_ 

9n 


rce 1 

Pr 


+  r  c(1-pt> 


(  4)  f  H-l )  0  f  1 


( 4>+n )  0  * 
2 


2 . 2-6 


+ 


30 
3  x 


e 


M 

3x 


The  boundary  conditions  are;  at  r\  =  0: 


<()w  =  0  (no  blowing  or  suction) 

(+)»  =  f'-l  (for  no  slip  f'  =  0) 

yw  w  w 
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<J>^  =  unknown  to  be  solved  for  2.2-7 

0w  "  8w-1  (sPecified) 

0^  =  8i 


as  ri  -+• 


<J>'  0 


2.2-8 


9+0 


Equations  2.2-5  and  2.2-6  are  coupled  and  therefore 
must  be  solved  simultaneously.  The  method  employed  here  is 
to  assume  C  =  1  for  the  first  iteration  which  uncouples  the 
momentum  equation  from  the  energy  equation.  The  momentum 
equation  can  be  solved  and  the  results  used  to  solve  the 
energy  equation.  The  ratio  C  can  now  be  expressed  as  a  func¬ 
tion  to  enthalpy  for  all  grid  points  using  the  values  of 
enthalpy  found  from  solving  the  energy  equation.  The  momen¬ 
tum  equation  is  solved  again  with  C  a  variable  and  then  the 
solution  to  the  energy  equation  is  obtained.  This  calculation 
is  repeated  until  the  results  converge. 

The  equations  are  solved  by  replacing  the  partial  deriva 
tives  with  respect  to  x  by  finite  difference  quotients  while 
retaining  the  derivatives  in  the  n  direction.  This  allows 
the  partial  derivatives  to  be  approximated  by  ordinary  diffe¬ 
rential  equations.  The  equations  are  solved  as  an  initial 


. 


value  problem. 


2.2.1  Solution  of  the  Momentum  Equation  The  first 
step  in  solving  the  momentum  equation  as  an  initial  value 
problem  is  to  determine 


cr  =  r  >dn  +  cw<*,w 

0 


The  next  step  is  to  solve  for, 


2.2-9 


=  W-i 

C 


2.2-10 


<f>*  =  /n  <£"  dn  + 

o 


2.2-11 


4> 


=  /n  cj)T  dn  +  (Pw 
0 


2 .2-12 


An  initial  guess  is  made  for  the  value  of  .  Since 

the  boundary  conditions  at  the  outer  edge  of  the  boundary 
layer  are  known  it  is  possible  to  vary  <p”  until  the  correct 
outer  boundary  conditions  are  me  1 1  The  solution,  however, 
can  be  obtained  more  rapidly  by  using  an  interpolation  pro- 

t 

cedure .  This  involves  obtaining  three  solutions  that  extend 
out  to  n„  within  the  bounds  of  4>'(n»)  +  1-0;  the  three  solu- 


.sraXdoiq  eu, 


ism  oifi  enotaibnoo  yifibnuod 
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tions  must  not  all  be  on  the  same  side  of  the  zero  axis. 

For  example,  two  solutions  can  have  4>  *  ( )  between  zero  and 
1.0  and  the  third  between  zero  and  -1.0.  Lagrangian  three 
point  interpolation  is  used  to  find  the  solution  which  meets 
the  outer  boundary  conditions  4>  *  ( Hoo )  =  0. 

Denoting  the  three  solutions  as 


<J>1 ,  4>2  >  $ 3 

4>  2»  $3 

4);’,  <f>’\  <f>” 


the  interpolated  solution  is  given  by 

<j>(n)  *  A1(f>1(n)  +  A24>2(n)  +  A3<f>3(n) 

<}> *  ( n )  *  A^<j)|(ri)  +  a202 (n)  +  a34>3 (n)  2.2-13 

cf)”  (n)  =  A^Cn)  +  A2 (n )  +  A3<t>^(n) 


where  the  coefficients  are 


A 


1 


4>^<n00) 

M>ilO  -  J  -  Cjivi 


' 


. 

(n)  A  *  (fr),  +  (pi  ;  ■  Oi)  » 
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a2 


*1(0 


■  [*J(n00)  -  <t>i (n^,)  ]  t*i(ri00)  -  *^ (!!„)] 


2.2-14 


♦{(loo)  $2  (loo) 


*  [*3(nm)  -  *i(n„)]  1 4>3 (n^)  -  ^(n.)] 


2.2.2  Solution  of  the  Energy  Equation  In  equation 
2.2-6,  6  and  0’  are  the  only  unknowns.  Since  the  equation 
is  treated  as  linear  it  is  possible  to  take  two  solutions  of 
equation  2.2-6  and  combine  them  to  meet  the  outer  boundary 
coundi t ion . 


For  simplicity  introduce  A  where 


2.2-15 


Hence  equation  2.2-6  becomes 


9(A) 

an 


(<Hn) 

2 


2.2-16 


Integrating , 


A  =  /n  A'  dn  +  Aw 
0 


2.2-17 


' 
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From  equation  2.2-15 


[A  -  si  c(1"p7>  U’+D*"] 


2.2-18 


and , 


2.2-19 


To  start  the  iteration  assume  some  value  of  0^.  In  the 
present  solution  a  value  of  0.2  was  used  and  the  solution  is 
denoted  as  0^(ri).  If  ©i(rioo)  is  greater  than  zero  a  smaller 
value  of  0^  is  tried;  if  it  is  less  a  larger  value  is  used. 
The  second  solution  is  denoted  as  ©2(0). 

The  two  solutions  are  linearly  independent  and  can 
be  added  to  give  the  general  solution  which  is  made  to  meet 
the  boundary  conditions. 


The  general  solution  is 


0(n)  =  A01(n)  +  b02(d) 


2.2-20 


The  boundary  conditions  at  the  outer  edge  are 


0  (rioo)  =  A01(noo)  +  B 0 2  (Ooo)  “  0 


e'(n„)  =  Ae^nJ  +  Be^nJ  =  o 


. 


28 


At  the  wall  the  boundary  conditions  are 


9w  “  A0Wl  +  B6 


w2 


However , 


0 


w 


Therefore , 

A  +  B  =  1 


and 


A 


-  02(ro 

-  e2<n„) 


The  correction  solution  is, 


0(n)  =  AG1(n)  +  ( l-A) e  2 ( n ) 


0' (n)  =  A0’ (n)  +  (1-A)0’ (n) 


2.2-21 


2 . 3  Boundary  Conditions 

The  velocity  and  temperature  profiles  are  obtained  from 
the  solution  of  the  boundary  layer  equations  with  appropriate 
boundary  conditions.  Two  sets  of  boundary  conditions  have 
been  used  in  this  thesis;  experimentally  determined  boundary 
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conditions,  and  the  classical  boundary  conditions  which  are 
defined  in  Section  2.3.2  below. 

Boundary  conditions  obtained  from  a  solution  of  the 
linearized  Boltzmann  equation  will  be  used  to  compare  with 
the  other  slip  resultSo 


2,3,1  Experimental  Results  Velocity  and  temperature 


profiles  can  be  calculated  using  experimental  values  of  tem¬ 
perature  and  velocity  at  the  wall  as  input  data,  together 
with  free  stream  static  pressure,  temperature,  and  velocity. 

2.3.2  Classical  Boundary  Conditions  A  semi  macro¬ 
scopic  argument  leads  to  the  following  relations  for  slip  and 
temperature  jump  [22,  23], 


U  (  x  ,  0  ) 


2  .  3-1 


and 


2 . 3-2 


where 


and 
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a0  =  -  — : 

2  a  8 


-  2-a  15 


2 


The  second  term  in  2.3-1  is  small  and  will  be  neglected 
in  this  analysis. 

Therefore , 


U  ( x  ,  0 )  =  H(ai  /rt  M) 
P  1  3  v 


2  „  3-3 


Equation  2.3-2  is  in  the  form  first  developed  by  Poisson 
and  equation  2.3-3  is  Maxwell's  equation  for  slip.  These 
boundary  conditions  will  be  referred  to  as  classical  boundary 
conditions  in  this  thesis. 

Slip  velocities  were  calculated  for  values  of  the  tan¬ 
gential  accommodation,  a,  equal  to  0.6,  0.7,  and  0.9  with 
the  thermal  accommodation  coefficient,  a,  taken  to  be  1.0. 

The  assumption  that  O  =  0.9  in  equation  2.3-3  corresponds  to 
Welander's  value  for  the  slip  coefficient  [24], 


2.3.3  Linearized  Bol t zmann  Solution  When  Rayleigh 


investigated  the  unsteady  flow  of  an  infinite  flat  plate  in 
an  incompressible  fluid,  he  was  able  to  draw  some  conclusions 
about  the  two-dimensional  flow  around  an  infinite  plate  by 
introducing  the  transformation  x  =  KUt,  where  U  is  the  plate 
velocity,  K  is  a  constant  taken  to  be  close  to  one  [25],  and 
x  is  the  distance  from  the  leading  edge. 


A  solution  of  the  linearized  time  variant  Boltzmann 
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equation  using  the  B.G.K,  model  for  collisions  has  been  ob¬ 
tained  by  C.  Cercignani  and  F.  Sernagiotto  [26],  It  is 
possible  to  transform  this  solution  to  a  steady  solution 
using  the  Rayleigh  transformation  and  thereby  obtain  a  slip 
boundary  condition.  It  must  be  noted,  however,  that  applying 
this  transformation  to  compressible  flow  may  give  only  quali¬ 
tative  results  [27],  Nevertheless  the  results  should  be  use¬ 
ful  to  compare  with  the  classical  boundary  conditions  which 
are  also  derived  from  a  somewhat  simplified  model. 

2.3.4  Incorporation  of  the  Boundary  Conditions  in  the 
Problem  The  technique  used  to  solve  the  boundary  layer  equa¬ 
tions  varies  slightly  depending  on  the  form  of  the  boundary 
conditions.  One  set  of  velocity  profiles  has  been  obtained 
using  experimental  velocity  and  temperature  boundary  condi¬ 
tions.  The  temperature  and  velocity  for  stations  along  the 
plate  are  read  into  the  program  as  starting  data  and  the  solu¬ 
tion  obtained  as  previously  described.  The  only  expression 

which  must  be  iterated  is  the  term  C  =  which  is  taken  to 

Poo  hoc 

be  equal  to  one  for  the  first  set  of  calculations.  In  sub¬ 
sequent  calculations  the  previous  calculation  is  used  to 
obtain  new  values  of  C.  It  was  found  that  after  two  itera¬ 
tions  of  C  the  velocity  profiles  remained  unchanged  for  the 
first  six  decimal  places. 

To  solve  the  boundary  layer  equations  with  the  classical 
boundary  conditions,  profiles  are  first  obtained  for  the  case 
of  zero  slip,  constant  wall  temperature,  and  C  =  1. 
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velocity  and  temperature  profiles  at  the  wall  are  used  to 
calculate  slip  and  temperature  jump  along  the  plate.  New 
values  for  C  are  obtained  as  in  the  previous  case.  This 
process  is  repeated  until  the  values  of  slip  and  tempera¬ 
ture  jump  have  converged  to  some  final  value.  In  practice, 
however,  the  number  of  iterations  is  limited  by  the  amount 
of  computer  time  available,,  It  was  found  that  five  iterations 
could  be  obtained  in  15  minutes  on  the  IBM  360-65,  and  that 
this  gave  results  which  converged  to  the  first  three  decimal 
places.  This  is  consistent  with  the  accuracy  of  the  inter¬ 
polation  procedure  (see  error  analysis  in  Chapter  V). 

2 . 4  Theory  of  the  Equilibrium  Temperature  Probe 

A  thermocouple  in  the  form  of  a  thin  wire  mounted  in  a 
free  molecular  flow  will  assume  an  equilibrium  temperature 
which  can  be  related  to  the  local  gas  temperature  in  terms 
of  speed  ratio,  S.  The  speed  ratio  is  obtained  directly  from 
impact  and  static  probe  measurements. 

The  relationship  between  equilibrium  wire  temperature, 
Tw,  and  local  gas  temperature,  T,  is  obtained  from  the  solu¬ 
tion  for  heat  transfer  between  a  gas  and  a  transverse  cylinder 
in  free  molecular  flow  which  has  been  derived  in  reference 
[28].  The  results  will  be  stated  here  and  the  details  of  the 
derivation  shown  in  Appendix  A. 

The  relationship  between  local  temperature,  wire  tem¬ 
perature,  and  speed  ratio  is, 
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( Z  x+Z  2 )  -  [Z1(S2  +  3)  +  Z2(S2+|-)] 


2.4-1 


+  2^/2 
PVma 


[ea(Tc4-Tt4)  -  Q] 


0 


where 


and 


Z1  =  TTe 


s2 


7TS  2  e 


7  S 2  S 2 

2  [  Iq  (~2~)  +  Il(V)] 


Q  =  internal  energy  input  to  the  cylinder  per  unit  area 
Tt  =  effective  temperature  of  the  surroundings 
Tc  =  wire  temperature 
T  =  local  temperature 
S  =  speed  ratio 

Iq  and  1-^  are  modified  Bessel  Functions  of  the  first 
and  second  kind.  Equation  2.4-1  is  simplified  further  since 
Q  =  0. 

A  calculation  in  Appendix  A  shows  that  the  radiation 
term  is  small  compared  to  the  other  terms. 

Neglecting  radiation  effects  the  equation  becomes, 
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[Z1(S2+3)  +  Z2 (S2+j)  ]  =  0  2.4-2 


CHAPTER  III 


APPARATUS 


3 . 1  The  Low  Density  Wind  Tunnel 

The  schematic  layout  and  a  photograph  of  the  low  density 
wind  tunnel  is  shown  in  Figures  1  and  2.  A  picture  of  the 
test  section  containing  the  flat  plate  model  is  shown  in 
Figures  3  and  4. 

The  tunnel  shell  is  24  inches  in  diameter  by  six  feet 
long  and  is  constructed  from  1/8  inch  stainless  steel.  The 
settling  chamber  is  a  12  inch  diameter  by  17  inches  long 
cylinder  fitted  with  a  removeable  nozzle  at  the  entrance  sec¬ 
tion.  Gas  enters  the  settling  chamber  through  a  length  of 
copper  line  with  a  1/4  inch  needle  valve  to  control  the  gas 
flow. 

The  system  is  roughed  with  two  Edwards  ES3000  mechani¬ 
cal  pumps  each  capable  of  handling  up  to  3000  cubic  feet  of 
gas  per  minute. 

Two  Edwards  B18  vapour  pumps  which  together  pump  5000 
liters  per  second  maintain  the  tunnel  at  test  conditions. 

The  pumps  are  oil  charged  and  have  a  water  cooling  system; 
a  12  inch  butterfly  valve  directly  above  each  diffusion  pump 
controls  the  pumping  speed. 

In  order  to  pump  still  larger  amounts  of  gas  a  22  inch 
diameter  hollow  cylinder  has  been  placed  downstream  from  the 
test  section  which,  when  filled  with  a  suitably  cold  fluid, 


. 
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removes  large  quantities  of  gas  by  cryogenic  pumping.  With 
a  tank  this  size  it  is  possible  to  attain  pumping  speeds  of 
the  order  of  several  hundred  thousand  liters  per  second 
[29].  A  convenient  arrangement  is  to  use  C02  as  the  working 
gas  and  maintain  a  cold  surface  with  liquid  nitrogen.  The 
liquid  nitrogen  level  may  be  controlled  by  an  automatic 
level  control  system  filling  from  a  large  storage  Dewar. 

Instrumentation  requiring  movement  in  the  axial  direc¬ 
tion  of  the  tunnel,  such  as  the  impact  probe,  can  be  mounted 
from  a  moveable  flat  plate  which  is  part  of  the  electron  gun 
housing  mechanism.  The  plate  is  propelled  from  outside  the 
tunnel  by  an  electric  motor  operating  through  a  mechanical 
feed  through  with  a  series  of  gear  reductions  (see  Figure  1) . 

3.1.1  The  Wind  Tunnel  Nozzle  A  schematic  drawing 
of  the  wind  tunnel  nozzle  is  shown  in  Figure  6.  It  was 
formed  out  of  two  pieces  of  brass  with  the  throat  area 
milled  from  a  solid  block  and  joined  to  a  length  of  brass 
pipe  which  forms  the  exit  section  of  the  nozzle. 

The  size  of  the  nozzle  is  governed  by  the  pumping  capa¬ 
city  of  the  Booster  Pumps.  Data  supplied  by  the  manufacturer 
indicates  a  capacity  of  2500  liters/sec.  at  a  pressure  of 
100  microns  for  the  two  B-18  pumps.  Making  an  allowance  for 
entry  losses  the  pumping  speed  was  assumed  to  be  2000  liters/ 
sec.  giving  a  nozzle  diameter  of  five  inches  for  flow  mach 


number  of  0.50. 


H„;b  3i  .  Id:'  '  ,:i£  ississl-?  '  •’ 


37 


3 . 2  Instrumentation 

3.2.1  The  Impact  Probe  A  view  of  the  impact  probe 
is  shown  in  Figure  5.  The  probe  was  made  of  stainless  steel 
in  the  shape  of  a  wedge  2  3/4  inches  long  by  1  1/8  inches 
wide.  The  opening  at  the  tip  is  1.0  inches  by  0.004  inches. 
This  allows  the  flow  to  be  sampled  to  within  0.002  inches 
of  the  model  surface. 

The  wedge  is  connected  by  a  four  inch  length  of  1/8  inch 
inside  diameter  tubing  to  a  temperature  compensated  NRC  801 
thermocouple  gauge  (an  analysis  of  the  temperature  dependence 
of  thermocouple  pressure  gauges  is  given  in  Appendix  B).  A 
digital  voltmeter  was  used  to  measure  the  output  of  the 
thermocouple  gauge,  giving  a  sensitivity  of  30  microvolts 
per  millitorr  in  the  pressure  range  of  50  to  100  microns. 
Repeatability  of  this  instrument  was  found  to  be  excellent. 

The  thermocouple  gauge  was  calibrated  against  a  Gran¬ 
ville  Phillips  capacitance  manometer  which  has  a  sensitivity 
of  10”“*  torr  if  the  voltage  response  is  measured  by  a  digital 
voltmeter . 

The  probe  was  traversed  through  the  boundary  layer  on 
a  mechanism  operated  from  outside  the  tunnel  by  means  of  a 
rotary  feed  through;  the  photograph  in  Figure  1  shows  the 
traversing  mechanism. 

A  dial  indicator  was  used  to  determined  the  vertical 
position  of  the  probe  relative  to  the  plate. 
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3.2.2  The  Static  Probe  Static  pressure  was  measured 
at  a  1/8  inch  diameter  static  tap  located  in  the  nozzle  wall 
1/2  inch  from  the  exit  (see  Figure  3).  An  NRC  801  thermo¬ 
couple  gauge  located  outside  the  tunnel  and  connected  to  the 
static  tap  by  a  short  length  of  stainless  steel  tubing  was 
used  to  measure  static  pressure.  The  digital  voltmeter  was 
also  used  to  measure  the  output  of  this  gauge  which  was  cali¬ 
brated  against  the  Granville  Phillips  Capacitance  manometer 

at  the  same  time  as  the  impact  pressure  gauge.  The  calibration 
curves  are  included  in  Figure  7. 

3.2.3  The  Equilibrium  Temperature  Probe  The  equili¬ 
brium  temperature  probe  is  shown  in  Figure  3.  It  consists 
of  a  0.003  inch  iron-cons tantan  thermocouple  wire  14  inches 
long,  supported  by  an  aluminum  arm.  The  ends  of  the  arm  are 
made  of  teflon  so  that  the  thermocouple  wire  is  insulated 
from  the  metal.  The  thermocouple  wire  is  held  in  tension 

by  passing  it  through  an  insulated  weak  spring  mounted  on  one 
of  the  arms  (Figure  3).  Iron-cons tantan  extension  wires  were 
soldered  to  the  ends  of  the  fine  wire  and  passed  through  the 
tunnel  wall  to  a  Leeds  and  Northrup  millivolt  potentiometer. 

The  thermocouple  junction  was  obtained  by  spot  welding 
a  pair  of  crossed  wires  and  then  trimming  the  ends  with  a 
sharp  razor  blade  under  a  microscope.  The  weld  is  about 
0.004  inches  in  diameter  which  gives  a  spatial  resolution 
compatible  with  that  of  the  impact  probe. 

The  Knudsen  number  for  the  equilibrium  temperature  probe 
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is  about  5.3  at  a  pressure  of  65  microns.  A  Knudsen  number 
greater  than  five  corresponds  to  the  free  molecular  flow 
regime  [ 30 ] . 

The  temperature  probe  was  mounted  in  the  support  pro¬ 
vided  for  the  impact  probe  and  its  position  measured  with 
the  dial  indicator  (see  Figure  4). 

The  plate  wall  temperature  was  measured  by  a  thermo¬ 
couple  attached  to  the  top  of  the  plate,  and  stagnation  tem¬ 
perature  was  measured  by  a  thermocouple  mounted  in  the 
stagnation  chamber.  Both  of  these  thermocouples  were  made 
from  0.022  inch  iron-cons  tan tan  wire.  All  three  thermocouple 
probes  were  read  on  the  same  potentiometer  using  a  Honeywell 
thermocouple  switch  and  an  ice  bath  for  a  reference  junction. 

3.2.4  The  Flat  Plate  Model  A  schematic  diagram  of 
the  flat  plate  model  is  shown  in  Figure  8.  Figure  3  and  4 
show  the  model  in  the  wind  tunnel  test  section. 

The  plate  was  made  of  aluminum  in  the  form  of  a  sharp 
wedge.  Two  small  compartments  were  milled  out  of  the  top  of 
the  wedge  and  two  thin  plates  welded  over  the  compartments 
in  such  a  way  that  the  compartments  are  flush  with  the  top 
of  the  plate.  Inlet  and  vent  lines  were  provided  for  filling 
the  plate  with  liquid  nitrogen;  a  hole  between  compartments 
keeps  the  liquid  level  the  same  on  both  sides.  Copper  lines 
are  used  to  carry  the  liquid  nitrogen  to  the  plate  from  out¬ 
side  the  tunnel  (Figure  4)  and  to  take  away  the  cold  gas. 

The  bottom  plate  surface  was  polished  smooth  and  the 
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leading  edge  brought  to  a  sharp  tip.  The  plate  was  placed 
in  the  tunnel  test  section  with  the  bottom  surface  parallel 
to  the  axis  of  the  nozzle  and  one  inch  above  the  nozzle 
centerline.  Positioning  of  the  plate  was  determined  from 
previously  conducted  nozzle  velocity  profile  surveys  (Figures 
9  and  10) . 

All  tests  were  conducted  on  the  bottom  side  of  the 
plate  as  this  insured  liquid  nitrogen  contact  with  the  test 
surface  for  the  cold  plate  tests  and  hence  less  likelihood 
of  temperature  fluctuations. 


CHAPTER  IV 


PROCEDURE 


4 . 1  General  Experimental  Procedure 

Experimental  investigations  were  carried  out  at  P  *  30 
microns  and  P  *  6 5  microns.  Two  series  of  tests  were  made 

at  65  microns,  one  with  an  uncooled  plate,  and  another  series 
with  the  plate  cooled  to  near  78°K  (liquid  nitrogen  tempera¬ 
ture)  . 

Impact  probe  measurements  were  taken  at  five  stations 
along  the  plate  and  up  to  1.0  inch  from  the  plate  surface. 

The  stations  along  the  plateare:  x  =  0.2,  0.4,  0.6,  0.8,  and 
1.0  inches,  measured  from  the  leading  edge.  Probe  measure¬ 
ments  were  obtained  next  at  11  wall  positions  starting  at 
the  leading  edge,  x  =  0,  and  moving  down  the  plate  in  0.1 
inch  intervals.  These  results  were  used  to  obtain  the  slip 
profile . 

The  equilibrium  temperature  probe  was  mounted  in  the 
support  used  for  the  impact  probe  measurements,  and  Tc 
was  obtained  at  the  above  mentioned  stations. 

Air  taken  into  the  tunnel  was  passed  through  a  desicant 
to  remove  water  vapour  which  would  have  formed  a  frost  on 
the  cold  plate.  There  was  no  visible  build  up  of  frost  on 
the  plate  for  the  duration  of  any  of  the  tests. 

The  cold  plate  temperature  was  maintained  constant  by 
filling  with  liquid  nitrogen  at  frequent  intervals  (about  10 
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minutes).  The  temperature  was  monitored  by  a  thermocouple 


mounted  on  the  plate. 

Since  most  of  the  tests  required  up  to  16  hours  to  com 
plete  some  fluctuations  in  conditions  were  encountered, 
especially  in  the  cold  plate  test  as  the  liquid  nitrogen 
passing  through  the  tunnel  wall  cooled  the  system  over  a 
period  of  hours.  In  order  to  reduce  fluctuation  errors  most 
of  the  experiments  were  run  at  least  twice  and  the  results 
averaged . 

4 . 2  Reduction  of  Experimental  Data 

If  the  flow  in  front  of  an  impact  probe  is  brought 
isen tropically  to  rest  the  relation  between  Pt,  P,  and  M  is. 


Pt  *  P(1  +  ’1~  M2)’^1 


4.2-1 


The  speed  ratio,  S,  defined  in  terms  of  mach  number  is 


S 


4.2-2 


Substituting  4.1-2  into  4.1-1  and  re-arranging 
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Speed  ratio  was  determined  from  the  measured  values 
of  P  and  Pt.  Local  temperature  was  obtained  from  equation 
2.5-2  using  the  speed  ratio  calculated  from  the  pressure 
measurements  together  with  the  wire  temperature.  Velocity 
can  then  be  found  from  the  definition  of  speed  ratio, 


U  =  /  2RT  S  4.2-4 


The  calibration  curve  in  Figure  7  was  used  to  reduce 
Pt  and  P  to  microns  pressure.  Tc  was  reduced  to  °K  using 

the  tables  in  reference  [31]. 

A  simple  computer  program  was  set  up  to  reduce  all  the 
experimental  data. 

4.2.1  Impact  Probe  Calibration  Pt  was  obtained  using 
a  probe  with  a  rectangular  orifice,  0.004  inches  high  and 
1.0  inches  wide.  It  was  felt  that  m  order  to  obtain  mean¬ 
ingful  results  good  spatial  resolution  would  be  required, 
especially  close  to  the  wall  where  gradients  are  large.  Cal¬ 
culations  show,  however,  that  the  response  time  for  a  circular 
open  ended  probe  of  say  0.004  inches  diameter  is  far  too 
long  to  make  measurements  practical.  Since  theoretical  cal¬ 
culations  of  response  time  in  the  transition  regime  may  be 
inexact,  small  impact  probes  varying  in  diameter  from  0.004 
inches  to  almost  .01  inches  were  made  and  tested  in  the  wind 
tunnel  at  P  *  70  microns.  The  response  time  of  the  smallest 

tube  was  of  the  order  of  1/2  hour  which  was  considered  im- 
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practical  for  the  large  number  of  points  to  be  measured. 

In  order  to  facilitate  the  reduction  of  the  data  it  is 
useful  to  keep  the  probe  in  either  free  molecular  flow  or  in 
continuum  flow.  At  P  =  65  microns,  X  *  0.03  inches,  which, 
if  the  crude  criteria  that  the  characteristic  dimension  must 
be  1/10  of  the  gas  mean  free  path  for  the  free  molecular 
probe  is  used,  would  require  a  circular  probe  of  0.003  inches 
diameter.  A  probe  of  this  dimension  would  be  impractical 
because  of  long  response  time.  Another  reason  for  rejecting 
a  free  molecular  probe  is  that  the  theory  used  to  reduce  the 
data  may  not  describe  the  flow  adequately  (see  reference  [32]). 
For  these  reasons  it  was  decided  to  use  a  continuum  probe  for 
the  measurements. 

In  order  to  insure  fast  response  it  was  decided  to  use 
a  probe  one  inch  wide.  The  problem  then  is  to  determine  for 
what  pressure  range  the  probe  is  in  the  continuum  regime 
since  it  is  not  known  what  dimension  should  be  taken  as 
characteristic  of  the  geometry.  A.  Sreekanth  [33]  has  shown, 
however,  that  perhaps  all  dimensions  of  the  probe  should  be 
considered  in  determining  whether  the  probe  is  in  the  con¬ 
tinuum  flow  regime  or  in  transition  flow.  For  this  reason  a 
test  was  conducted  to  find  the  continuum  limit  of  the  wedge. 

The  impact  pressure  for  static  pressures  ranging  from  90 
microns  to  17  microns  were  measured  and  compared  to  results 
obtained  with  a  circular  probe  with  a  0.187  inch  diameter 
orifice.  The  results  are  shown  in  Figure  11  with  impact 
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pressures  agreeing  very  well  down  to  a  static  pressure  of 
50  microns.  Below  this  pressure  the  wedge  pressures  are 
lower  than  the  circular  probe  pressures.  This  would  indi¬ 
cate  that  both  probes  are  in  the  continuum  down  to  50 
microns  where  it  appears  that  the  circular  probe  enters  the 
transition  regime.  At  some  lower  value  of  pressure  the 
wedge  probe  will  be  in  transition  flow. 


CHAPTER  V 


RESULTS  AND  DISCUSSION  OF  EXPERIMENTAL  AND 
THEORETICAL  INVESTIGATION 


5 . 1  Slip  and  Temperature  Jump  Boundary  Conditions 

5,1,1  Results 

5 . 1 o 1 . 1  Measurement  of  Slip  Velocity  on  a  Cold 
Plate  Slip  velocity  as  a  function  of  distance  from  the 
leading  edge  was  experimentally  measured,  These  results  are 
shown  in  Table  1  and  plotted  in  Figure  12 , 

Five  Theoretical  determinations  of  slip  velocity  are 
shown  in  Figure  12  for  comparison.  The  curves  marked  classi¬ 
cal  boundary  conditions  represent  the  end  result  of  the 
iterative  solution  of  the  boundary  layer  equations  using  the 
classical  slip  condition  and,  where  indicated,  the  tempera¬ 
ture  jump  boundary  condition  given  in  Section  2.4.2.  The 
boundary  layer  equations  have  been  solved  for  a  constant 
Prandtl  number  of  0,76.  Summarizing  the  theoretical  results, 
they  include: 

(a)  A  solution  of  the  classical  slip  condition  with  tempera¬ 
ture  jump,  and  with  O  =  0,9,  ot  =  1,0.* 

(b)  The  Rayleigh  solution  discussed  in  Section  2.4.2 


*  Unless  otherwise  indicated  the  thermal  accommodation  co¬ 
efficient  will  be  taken  to  be  1.0, 


(c)  A  solution  of  the  classical  slip  condition  for  no  tempera 
ture  jump,  with  o  =  0.9,  and  > 0  ^  =  0.5. 

U  oo 

(d)  A  solution  of  the  classical  slip  condition  with  no  tem¬ 
perature  jump;  0  =  0,9  and  — =  1.0. 

Uqo 

(e)  A  solution  with  classical  slip  and  temperature  jump; 

a  =  0.6  and  y.(ff-*-Q)  =  lo0. 

U  00 


5 . 1 c 1 , 2  Measurement  of  Slip  and  Temperature  Jump 
on  an  Uncooled  Plate  at  P  =  30.4  Microns  Figure  13  shows 
the  experimental  and  theoretical  results  for  slip  over  an 
uncooled  plate  at  P  =  30=4  microns,  and  Figure  14  shows  the 
experimental  and  theoretical  results  for  temperature  jump. 

The  experimental  data  for  slip  and  temperature  jump  is 
summarized  in  Table  2, 

The  theoretical  slip  profiles  obtained  from  a  solution 
of  the  boundary  layer  equations  for  a  constant  Prandtl  number 
of  0.71  include : 

(a)  A  solution  of  the  classical  slip  condition  with  tempera- 

A  n  j  U(0,0)  "|  A 

ture  jump;  O  =  0,9  and  - — 1 -  -  1  =  0, 

Uqo 

(b)  A  solution  of  the  classical  slip  condition  with  no  tem- 

„  „  ,  U  ( 0 , 0 )  _  A 

perature  jump;  O  -  0,9  and  ^  ~  1.0, 

(c)  A  solution  of  the  classical  slip  condition  with  tempera¬ 
ture  jump;  O  =  0.7  and  — ~  =  1*0 

Uqo 

(d)  The  Rayleigh  solution. 

The  theoretical  temperature  jump  curve  in  Figure  14  is 

the  curve  obtained  in  conjunction  with  the  solution  in  (a) , 
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5 . 1 . 1 . 3  Measurement  of  Slip  and  Temperature  Jump 
on  an  Uncooled  Plate  at  P  =  65.0  Microns  A  test  similar  to 
the  one  in  Section  5.1. 1.2  was  conducted  at  P  =  65.0  microns. 
Slip  profiles  are  shown  in  Figure  15  and  temperature  jump 
curves  in  Figure  16.  The  experimental  data  is  tabulated  in 
Table  3. 


The  theoretical  slip  profiles  include: 


(a)  A  classical  slip  profile  with  temperature  jump  effects; 


(b) 


(c) 


(d) 


U(Q>  =  1.0  and  a  =  0.6. 

U  r 


'  00 


A  classical  slip  profile  with  no  temperature  jump  effect; 
U ( 0  >  =  0.5  and  a  =  0.9. 


U 


00 


Classical  slip  with  no  temperature  jump  effect;  ?  ^  = 

U  r 


'00 


1.0  and  O  =  0.9. 

Classical  slip  with  temperature  jump  effects;  U  =  1.0 

^  t 


'00 


and  a  =  0.9. 


(e)  The  Rayleigh  solution. 

The  theoretical  temperature  jump  profile  in  Figure  16 
is  obtained  in  conjunction  with  the  solution  in  (d)  . 


5.1.2  Discussion 


5. 1.2.1  Experimental  Results  Examination  of  the 
slip  profiles  show  an  apparent  leveling  off  of  the  curves  at 
approximately  0.3  inches  in  the  cold  plate  test,  and  at  0.6 
to  0.7  inches  in  the  uncooled  plate  tests.  Values  of  slip 
velocity  after  1.0  inch  are  =  0.279  and  =  0.336  for 
the  uncooled  plate  at  30.4  microns  and  65.0  microns  respec- 
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tively,  and  — —  -  0.356  for  the  cooled  plate  at  65.0  microns 


00 


These  values  are  somewhat  higher  than  the  theory  pre¬ 
dicts.  Some  higher  order  terms  which  are  important  near 
the  leading  edge  have  been  neglected  in  the  theory  and 
inclusion  of  these  would  tend  to  increase  the  theoretical 
values  somewhat. 

Small  changes  in  static  pressure  would  have  a  large 
influence  on  the  experimental  slip  values,  particularly  when 
the  slip  velocities  are  small  as  can  be  seen  from  equation 
4.2-3.  For  example,  an  increase  of  static  pressure  of  0.5 
microns  in  the  run  at  30.5  microns,  and  1.6  microns  in  the 
run  at  65.0  microns,  is  enough  to  make  the  experiments  coincide 
with  theory.  Experimental  curves  have  been  calculated  with 
these  changes  in  static  pressure  and  are  shown  in  Figures 
13  and  15 . 

Static  pressure  was  assumed  to  be  constant  over  the  plate 
and  equal  to  the  static  pressure  measured  1/2  inch  from  the 
nozzle  exit.  An  increase  in  static  pressure  could  occur  due 
to  such  effects  as  spreading  of  the  free  jet  and  influence 
of  the  model  and  impact  probe  in  the  stream.  Any  one  of  these 
could  account  for  the  small  pressure  change  needed  to  account 
for  the  difference  between  experiment  and  theory.  It  is 
intfirGSting  to  note  that  in  the  uncooled  tests,  the  measured 
value  of  —  is  0.903  at  x  =  0  when  P  =  30.4  microns,  and 

Uoo 

JL  =  0  893  at  P  =  65.0  microns.  These  values  are  somewhat 
Uoo 

higher  than  some  theories  predict.  For  example  references 
[26]  [34]  and  [35]  obtain  =  0.5  at  the  leading  edge  when 
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a  =  1.0 

The  experimental  temperature  jump  profiles  level  off 
after  approximately  0.2  inches  from  the  leading  edge  and 
differ  with  theory  by  less  than  1%  at  x  =  1.0  inch.  The 
data  clearly  indicates  that  non  adiabatic  wall  conditions 
exist  in  the  test  region;  for  example,  Tw  =  301. 4°K  and 
T^  =  291. 8°K  for  the  P  =  30.4  micron  case,  and  Tw  =  300. 5°K 
and  T^  =  296. 5°K  at  P  =  65.0  microns. 

5. 1.2. 2  Theoretical  Results  Comparison  of  the 
theoretical  profiles  with  experiment  shows  that  theory  pre¬ 
dicts  lower  values  of  slip  and  that  the  best  correlation 
occurs  for  the  solutions  having  lower  values  of  accommodation. 

E.  Merlic  [12]  measured  the  tangential  accommodation 
of  air  on  an  aluminum  cylinder  and  found  that  a  =  0.9  on  a 
surface  left  in  a  chamber  at  0.01  microns  for  one  to  three 
days  before  carrying  out  tests.  He  found,  however,  that  if 
the  chamber  was  maintained  at  250  microns  for  six  days  before 
tests,  a  value  of  a  =  0.6  was  obtained.  This  was  thought  to 
be  due  to  the  adsorption  of  an  oil  film  on  the  aluminum  sur¬ 
face,  caused  by  back  migration  from  the  diffusion  pump.  The 
oil  film  is  more  "specular"  than  the  metal  surface  and  hence 
lowers  the  value  of  the  accommodation  coefficient.  This 
condition  of  the  surface  is  similar  to  that  of  the  aluminum 
plate  in  the  low  density  wind  tunnel;  the  plate  was  cleaned 
after  mounting  but  the  tunnel  was  pumped  for  a  few  days  before 
actual  tests  began.  Furthermore,  the  physical  size  of  the 
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test  section  makes  it  necessary  to  mount  the  model  near  the 
opening  to  one  of  the  diffusion  pumps,  increasing  the  possi¬ 
bility  of  oil  contamination.  Baffles  placed  on  the  opening 
to  the  diffusion  pumps  reduce  the  problem  but  do  not  com¬ 
pletely  eliminate  it. 

Oil  contamination  of  the  surface  was  increased  in  the 
cold  plate  tests  where  the  plate  acts  as  a  very  effective 
cold  trap.  After  the  cold  plate  tests  a  slight  oil  film 
could  be  detected  by  running  a  finger  over  the  plate  surface. 

Theoretical  solutions  of  classical  slip  with  no  tempera¬ 
ture  jump  have  been  obtained  with  starting  values  of  =  1.0 

^00 

and  0.5  at  the  leading  edge  to  observe  the  change  in  the  slip 
profiles  caused  by  varying  the  starting  values.  A  very 
important  observation  is  that  at  P  =  65.0  microns  the  uncooled 
plate  slip  velocities  are  almost  the  same  after  0.15  inches 
from  the  leading  edge  for  both  starting  conditions.  The  same 
effect  is  observed  in  the  cold  plate  profiles  with  both 
curves  merging  at  approximately  X  =  0.25  inches. 

Another  interesting  point  is  that  inclusion  of  the  tem¬ 
perature  jump  boundary  condition  does  not  have  a  very  signi¬ 
ficant  affect  on  the  slip  profiles.  As  expected,  however, 
the  addition  of  this  boundary  condition  becomes  more  important 
at  lower  pressures.  This  can  be  seen  by  comparing  Figure  14 

and  Figure  16 . 

The  Rayleigh  solution  included  in  the  theory  compares 
well  with  classical  theory  for  an  uncooled  plate  if  the 
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accommodation  coefficient,  a,  is  assumed  to  be  0.9.  Since 
the  Rayleigh  solution  starts  with  — -  =  0.5  at  the  leading 

^  OO 

edge  it  will  compare  best  with  the  classical  solution  having 

TT”  =  0*^  at  x  =  0  (see  Figure  15).  The  Rayleigh  solution 
^00 

does  not  compare  as  well  with  the  classical  theory  for  flow 
over  a  very  cold  plate  (see  Figure  12) . 

5 . 2  Velocity  and  Temperature  Profiles 

5.2.1  Results 

5 . 2 . 1 . 1  Velocity  and  Temperature  Profiles  on  a 
Cold  Plate  Velocity  and  temperature  profiles  for  flow  over 
a  cold  flat  plate  are  summarized  in  Tables  4-8  with  velocity 
versus  distance  from  the  plate  surface  at  five  different 
plate  positions  plotted  in  Figures  17  to  21.  Temperature 
profiles  are  shown  in  Figures  22  to  26.  Two  sets  of  theore¬ 
tical  curves,  one  using  as  starting  values  the  previously 
mentioned  experimental  slip  and  temperature  jump  profiles  with 
corrections  for  small  changes  in  P,  ,  and  T^,  and  one  using 
the  classical  boundary  conditions  (a  =  0.90)  are  included  for 
comparison  with  the  experimental  curves.  The  plate  wall  tem¬ 
perature  is  considered  considered  constant  at  82.7  K  and 

the  value  of  the  free  stream  Mach  number  is  0.63. 

5.2. 1.2  Velocity  and  Temperature  Profiles  on  an 
Uncooled  Plate  Results  for  experiments  at  P  =  29.0  microns 
and  P  =  64.0  microns  are  tabulated  in  Tables  9-13  and  14-18 


. 
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respectively.  Velocity  profiles  are  shown  in  Figures  27  to 

30  and  Figures  35  to  39,  and  temperature  profiles  in  Figures 

31  to  34,  and  Figures  40  to  44.  Two  theoretical  determina¬ 
tions  using  the  boundary  conditions  mentioned  in  5. 2. 1.1 
have  been  included. 

5.2.2  Discuss  ion  A  general  observation  of  velocity 
and  temperature  profiles  indicates  a  large  discrepancy  between 
theory  and  experiment  at  x  =  0.2  inches  with  considerable 
improvement  in  the  correlation  after  1.0  inches  in  the  un¬ 
cooled  plate  tests.  At  1.0  inches  Re  *  11  for  P  =  30.4 

X 

microns  and  Rex  ~  24  for  P  =  65.0  microns.  Good  agreement 
would  not  be  expected  in  the  region  near  the  leading  edge 
considering  that  the  order  of  magnitude  analysis  on  the 
governing  theoretical  equations  in  Chapter  II  assumed 
ReL  »  0(1). 

5 . 3  Error  Analysis 

5.3.1  Experimental  Accuracy 

5 . 3 . 1 . 1  Resolution  of  the  Static  and  Impact  Probes 
The  instrument  used  to  measure  impact  and  static  pressures 
has  a  resolution  of  0.3  microns  which  results  in  an  error  of 
+  0.3%  to  +  1%  in  the  pressure  data.  This  error  has  been 
reduced  by  taking  more  than  one  pressure  reading  and  averag¬ 
ing  the  results.  Since  the  error  caused  by  the  resolution 
of  the  instrument  will  be  larger  at  lower  pressures,  the  error 
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in  equation  4.2-3  has  been  calculated  for  P  %  30  microns  and 
P^_  »  P.  The  static  pressures  have  been  arrived  at  by  averag¬ 
ing  up  to  10  pressure  readings  and  if  they  can  be  considered 
correct,  the  maximum  error  in  equation  4.2-3  is  6.5%.  The 
magnitude  of  the  error,  however,  will  be  much  less  than 
this  for  most  of  the  tests. 

5 . 3 . 1 . 2  Accuracy  of  Temperature  Measurements 

The  accuracy  of  the  equilibrium  wire  temperature  measurements 
vary  from  0.33°C  to  0.20°C  depending  on  the  magnitude  of  the 
temperature.  For  example,  at  temperature  of  the  order  -109°C 
the  tables  in  reference  [31]  can  be  read  to  0.33°C.  This  is 
considered  to  be  acceptable  when  taking  into  account  other 
approximations . 

Values  of  local  temperature  are  obtained  from  expression 

2.5- 2.  A  calculation  shows  that  the  effect  of  radiation  is 
negligible  and  therefore  it  has  not  been  considered  in  the 

calculations  (see  Appendix  A) . 

The  maximum  error  in  speed  ratio  was  calculated  in  the 
previous  section  to  be  6.5%.  Using  typical  values  of  Tc  = 

301 °K  and  S  -  0.10  the  error  in  T  calculated  from  equation 

2.5- 2  is  0.07%. 

5. 3. 1.3  Error  in  Reading  the  Dial  Indicator  Read¬ 
ings  were  taken  at  various  positions  in  the  boundary  layer  with 
a  dial  indicator  capable  of  being  read  to  0.001  inches.  How¬ 
ever,  errors  could  result  from  a  certain  amount  of  back  lash 

of  the  traversing  mechanism  (see  Chapter  III). 
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5. 3. 1.4  Other  Possible  Sources  of  Error  Other 
possible  sources  of  error  not  previously  mentioned  are  the 
spreading  of  the  nozzle  jet  and  the  influence  of  the  model 
on  the  flow.  It  is  necessary  to  treat  these  effects  simul¬ 
taneous  ly  . 

S.  Gordon  [32]  investigated  the  flow  of  gas  over  a  two 
inch  flat  plate  lying  one  inch  in  front  of  the  nozzle  exit. 

He  observed  that  the  plate  influenced  the  flow  up  to  approxi¬ 
mately  10  mean  free  paths  in  front  of  the  leading  edge.  To 
complicate  matters  he  also  noted  that  the  centerline  speed 
ratio  had  dropped  by  12%  three  inches  from  the  nozzle  exit 
which  he  attributed  to  the  spreading  of  the  nozzle  jet.  From 
Gordon's  observations  it  would  seem  necessary  to  place  the 
model  far  enough  from  the  nozzle  exit  so  as  not  to  disturb 
the  flow  in  the  nozzle,  and  at  the  same  time,  close  enough  to 
minimize  the  spreading  of  the  nozzle  jet  on  the  area  of  the 
plate  tested.  For  the  present  tests  the  model  was  placed  0.22 
inches  from  the  nozzle  exit  which  corresponds  to  about  eight 
mean  free  paths  at  65.0  microns.  There  was  no  measureable 
drop  in  centerline  velocity  for  the  region  studied.  Undoub¬ 
tedly  the  presence  of  the  model  influenced  the  flow  in  front 
of  the  leading  edge  but  its  affect  has  not  been  estimated. 

5.3.2  Accuracy  of  the  Numerical  Analysis  The  errors 
involved  in  using  the  integration  and  finite  difference  for¬ 
mulas  are  summarized  in  Appendix  C.  The  accuracy  of  the 
numerical  technique  for  the  compressible  boundary  layer  with¬ 
out  slip  has  been  analyzed  in  reference  [22]. 
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The  solution  at  a  station  is  obtained  by  getting  three 
solutions  that  extend  all  the  way  to  ri^,  within  the  bounds  of 
4>f  (O  +  Kf  where  K'  has  been  taken  as  1.0  in  the  analysis. 

This  bound  was  checked  by  running  a  program  with  K'  =  0.95 
which  showed  a  negligible  change  in  the  results.  The  step 
size  in  the  x  direction  was  taken  to  be  0.02  for  the  first  0.1 
inches  and  then  increased  to  0.04  for  the  remaining  stations. 

The  error  in  the  results  are  not  necessarily  decreased  by  tak¬ 
ing  smaller  step  sizes  in  the  x  direction  as  the  important 
parameter  is  .  Oscillations  in  the  curves  can  be  prevented 
by  keeping  £  25  [21].  This  criterion  was  proven  to  be 

correct  as  it  was  found  that  oscillations  began  to  appear  for 

~  28. 

Ax 

The  change  in  r|  was  chosen  to  be  0.1  and  rioo  =  5.2. 

These  values  are  governed  largely  by  computer  storage  space; 
however,  it  was  found  that  an  increase  to  Ooo  =  5.7  only  affected 
the  fourth  or  fifth  decimal  place  in  the  velocity  results. 

Since  the  analysis  is  accurate  to  the  third  or  fourth  decimal 
place  this  difference  is  not  significant. 

The  Prandtl  number  was  taken  to  be  0.71  for  the  uncooled 
tests  and  0.76  for  the  cooled  test.  The  effect  of  Prandtl 
number  on  the  results  was  checked  by  changing  the  value  from 
0.76  to  0.73  in  the  cold  plate  theory.  This  caused  the 
boundary  layer  thickness  to  vary  by  0.0006  inches  at  x  = 

0.02  inches  and  0.0045  inches  at  x  =  1.0  inches. 
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CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 


A  theoretical  and  experimental  study  has  been  conducted 
in  the  transition  flow  regime.  The  conclusions  from  the 
investigation  are  listed  as  follows: 

(a)  The  numerical  analysis  of  the  compressible  boun¬ 
dary  layer  on  an  uncooled  plate  with  experimentally  deter¬ 
mined  boundary  conditions  indicates  that  the  boundary  layer 
equations  are  valid  for  x  1.0  inches  and  Kn  ~  0.11.  At 

x  =  1.0  inches  and  P  *  30  microns  the  value  of  the  local 
Reynolds  number  is  approximately  11. 

(b)  Theory  predicts  lower  values  of  slip  than  were 


experimentally  determined. 

(c)  Experimental  results  give  a  value  of 
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0.9  at 
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the  leading  edge,  contrary  to  some  theories  which  obtain 
—  =  0.5  for  diffuse  reflection.  The  starting  value  for 

Uqo  U00 

at  the  leading  edge  does  not  significantly  affect  the  slip 

profiles  after  approximately  three  mean  free  paths  from  the 


leading  edge. 

(d)  The  classical  slip  condition  as  determined  from  an 
iterative  solution  of  the  boundary  layer  equations  has  been 
compared  to  Cercignani  and  Sernagiotto's  time  varient  solution 
of  the  Boltzmann  equation,  transformed  to  a  two-dimensional 
steady  problem.  The  theories  agree  very  well  for  a  tangen¬ 
tial  accommodation  of  0.9  in  the  classical  theory. 
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(e)  For  Kn  ~  0.6  the  inclusion  of  the  temperature  jump 
boundary  condition  does  not  change  the  slip  velocity  profile 
appreciably;  for  Kn  ~  0.11  the  change  is  more  significant 
but  still  small. 


(f)  The  use  of  a  rectangular  shaped  impact  probe  gives 
better  resolution  and  can  measure  total  pressures  closer  to 
the  wall  than  a  circular  orifice  having  a  similar  response 
time  . 


(g)  A  low  density  wind  tunnel  has  been  designed  and 
ins trumen tated  which  enables  the  study  of  transition  flows. 
The  tunnel  operates  under  stable  conditions  with  very  little 
evidence  of  oil  back  migration. 
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THEORY  OF  THE  EQUILIBRIUM  TEMPERATURE  PROBE 
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The  basic  assumption  in  the  derivation  of  the  solution 
for  heat  transfer  between  a  gas  and  a  transverse  cylinder  in 
free  molecular  flow  is  that  the  impinging  molecules  have  a 
classical  Maxwellian  velocity  distribution  superimposed  on 
a  uniform  mass  velocity. 

For  a  monatomic  gas  the  translational  molecular  energy 
incident  on  an  element  of  area  dA  inclined  at  an  angle  6  with 
respect  to  the  stream  mass  velocity  U  is, 

dE  t  =  n(^p  +  ip  kT )  dA  A-l 


where  k  =  Boltzmann  constant,  and, 


l+i£T  s  sin  6  [1  +  erf  (s  sin  e)]es2sin2e 

ijj  =  1  +  - — — -  ■  '  2  !  2  n  ^ 

1+  /rT  s  sin  0[1  +  erf  (s  sin  9)]es  sin 


n  is  the  number  of  molecules  striking  a  unit  surface 
area  per  unit  time  and  can  be  expressed  by 


NVmX 
n  =  - 

2/tT 


A-  3 


The  energy 


incident  from  the  opposite  side  is  given  by 


' 
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dE 


t 


+  \\)'kT)dA 


A- 4 


wi  th 


and 


2/tt 


A-  5 


1  — s  sin  6 [ 1  -  erf (s  sin  6)]es^sin^® 

r-  i  +  — £ -  a- 6 

o  9 

1-  /T  s  sin  0  [  1  -  erf  (s  sin  9)]es  s^n  ® 


For  a  diatomic  gas  the  rotational  component  of  the 
molecular  internal  energy  must  be  added  to  the  translational 
component.  Considering  a  diatomic  molecule  whose  structure 
is  assumed  to  be  that  of  a  dumbbell 

Er  =  kT  A- 7 

The  total  molecular  energy  becomes 


dE 


i 


+  (l|j+l)kT]dA 


A-  8 


dE'  =  n  "  [-J-  +  (ip"+l)kT]dA 


A-  9 


• 

. 
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N  ow , 


A- 10 


where  a  is  defined  as  the  thermal  accommodation  coefficient 
and  Ew  is  the  rate  of  re-emitted  molecular  energy  carried  by 
the  scattered  stream  if  it  is  in  equilibrium  with  the  surface 
temperature  Tw . 

The  rate  of  re-emitted  energy  for  a  diatomic  gas  is 
given  by, 


dEw  =  3nkTwdA  A-ll 

Combining  A-10  and  A-ll  the  energy  transported  from  the 
front  of  the  element  by  re-emitted  molecules  becomes 

dE  =  (dE.  +  nkTdA)  (1-a)  +  3ankT  dA  A-12 

r  t  *- 

where  Tc  is  the  surface  temperature  of  the  cylinder. 

Similarly  the  energy  transported  from  the  rear  side  of 
dA  is  given  by 


dE'  =  (dE£  +  n'kTdA)  (1-a)  +  3an'kTcdA  A-13 

For  a  cylinder  of  length  L  and  radius  r  the  elemental 


area  is 


eemoosd  89IUosio«  V«'  **•  °  ’"* 


' 


dA 


rLdQ 


6  6 


The  energy  balance  is 


y  2  ±2 

2arL{/2  n[^-  +  (i|M-l)kT]d0  +  j2  n'[-~  +  (i|T+l)  kT  ]  d0  } 

0  0  2 

7T  TT 

2  2 

-  6arLkT C(J  nd0  +  /  n^d0)  A-14 

0  0 

+  2TTrL[Q  -  ea( Tc4  -  Tt4)  ]  =  0 

where  Q  is  the  internal  energy  input  to  the  cylinder  per  unit 
area  and  it  has  been  assumed  gray  body  radiation  occurs 
between  the  cylinder  and  the  surroundings.  Tt  is  the  effec¬ 
tive  temperature  of  the  surroundings. 

The  integrals  in  equation  A-14  have  been  solved  in 
[28],  The  solution  then  becomes 


3X(Z1+Z2} 


[Z1(s2+3)  +  Z2 (s2+j) ] 


A=  1 5 


2  TT 


3/2 


PVma 


[ea(T, 


-  Tt*> 


-  Q) ]  =  0 


8,1000  noisslba,  <bcd  <618  baauae*  **»<!  ■*«'  11  b“  ** 

' 

5 


where 


2 


s 


2 


2 


Z  ]_  =  7T  e 


and 


Z 


2 


Iq  and  1^  denote  modified  Bessel  functions  of  the  first 
kind  and  zero  and  first  order  respectively. 

[36]  gives  a  series  expansion  of  the  modified  Bessel 
function  for  small  arguments.  Both  series  may  be  truncated 
after  the  second  term  as  the  error  in  Iq  and  1^  introduced 
by  neglecting  higher  order  terms  is  0.36%  and  0.03%  respec¬ 
tively  for  a  speed  ratio  of  1.0.  Since  the  maximum  speed 
ratio  is  about  0.7  this  error  is  reduced  considerably  further. 

A.l  Radiation  Effects 


To  analyze  the  magnitude  of  the  radiation  term  in  A-15 


a  simple  calculation  is  carried  out. 


Ass  uming : 

S  =  0.70 

Tc  =  30 0 

Tt  =  -190° 

p  =  0.0056  lb/in2  (*  30  microns) 


•• 
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e  =0.3 

a  =  1 

Vm  =  238  f  t/sec 

substituting  these  values  into  A-15  gives 


T 

3-~-(  3.76)  -  10  c  06  +  0.5  =  0 


For  T c  =  30  °  C , 

T  =  34 . 1 °  C . 

Neglecting  radition  effects  gives 

T  =  34 . 0  °  C 

This  gives  an  error  of  .29%  which  is  less  than  the 
other  errors  involved  in  measuring  temperature. 

Therefore  neglecting  radiation  effects  the  final  equa¬ 
tion  becomes 


3^-(Z1+Z2)  -  [Z1(s2  +  3)  +  Z 2  ( s 2  +~ )  ]  =  0 


A- 16 
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APPENDIX  B 


PRESSURE  RESPONSE  TO  THERMOCOUPLE  WALL 


TEMPERATURE  VARIATION 
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The  thermocouple  gauge  works  on  the  principle  that  the 
heat  loss  to  a  hot  wire  in  free  molecular  flow  is  proportional 
to  the  pressure  in  the  thermocouple  head. 

Since  the  important  modes  of  heat  transfer  are  conduc¬ 
tion  by  gas  molecules  and  radiation,  the  energy  balance  will 
consist  of  the  addition  of  these  effects.  However,  radiation 
effects  only  become  important  at  low  pressures  [37];  for 
example,  at  P  8  0.01  torr  with  a  filament  temperature  of 
100°C  and  wall  temperature  at  20°C,  the  heat  transfer  due  to 
conduction  is  8.87  x  10-2  watts/cm2  and,  that  due  to  radiation 
is  6.80  x  10“^  watts/cm2.  As  the  tests  conducted  at  P  z  30 
microns  will  have  a  still  larger  relative  heat  transfer  by 
conduction,  radiation  will  be  neglected  in  this  simple  analy¬ 
sis. 

Therefore , 


Et  -  PK0a(|H-)1/2  (Tc  -  Tgas) 
gas 

where  Kq  is  the  free  molecular  conductivity  of  air  and 

Tgas  ”  Tw 

To  obtain  the  expression  for  response  due  to  temperature 
change,  assume  two  different  sets  of  values  of  Tw  and  P  which 


give  the  same  value  for  Et. 
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Equating  and  re-arranging, 


(Tp  -  Tw_) 


'  m  '  T  _  T 

iW2  c  wj 


Substituting  Tc  =  150°C,  TW2  =  20°C,  P2  =  30  microns, 

Tw^  =  10°C,  the  value  for  P-j_  is  27,4  microns,  This  represents 

an  error  of  8.66%  which  is  significant.  When  dealing  with 
the  cooled  plate,  a  temperature  change  of  this  magnitude  could 
easily  occur.  Because  of  the  difficulty  in  maintaining  con¬ 
stant  thermocouple  wall  temperature  it  was  deemed  necessary 
to  have  a  temperature  compensated  thermocouple  gauge  for  the 
experiments.  The  National  Research  Corporation  have  a  thermo¬ 
couple  gauge  which  is  temperature  compensated,  and  it  was 
decided  to  use  the  NRC  801  thermocouple  gauge. 

The  outline  of  a  design  of  a  sensitive  pirani  gauge  is 
given  in  reference  [39].  This  gauge  has  a  sensitivity  of  40 
millivolts  per  millitorr  and  should  give  better  resolution 
than  the  thermocouple  gauge  used.  However,  it  is  not  readily 
available  and  may  require  special  corrections  for  temperature 


variations . 
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ERROR  IN  NUMERICAL  EXPRESSIONS 
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The  errors  involved  in  the  various  experssions  in  the 
numerical  analysis  are  listed  here. 

C o 1  Finite  Difference  Representations 
Two  Points : 


~  ^M-l 
9x 


The  error  is , 


(Xm-xm-i)  a2(j)  (Q 


dx‘ 


where  E,  is  some  value  of  x  in  the  interval  X^  - 
Three  Points : 


3  x 


(XM  XM-1')  (XM"XM-2) 


]cf) 


M 


(XM-XM-2> 

(  XM_XM- 1  ^  ('XM-l"XM-2) 


^M-l 


_ (  xm~xM— l) _ 

+  (xm-%-2^xM-1"XM-2^ 


The  error  is , 


■  ’ 


7  4 


(XM  XM-1}  (XM  XM-2}  a34>(£) 

6  8x3 


C  o  2  Integration  Expressions 


C . 2 . 1  Extrapolation  Formulas 


(cr>N)E  =  (cr)N_x  +  ^[55(C 


-  59  (C4>”  )^_2  +  37(Cc})")n-3  "  9(C4)")^_4] 


where  An  is  the  step  size  and  subscript  E  denotes  extrapolation. 
The  error  is 


ex  <  +  (An)5  [C<t>" (5)  1 v 


The  formulas  for  <p '  and  cf)  are: 


*N-1 


+ 


An 

24 


[55^_i 


39^.2  +  37^_3  -  94>JJ_4] 


F  <  |  251 

E2  -  720 


(An)5  4>VI(5) 


and  , 
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^N)E  "  ^N-l  +  An^N-l 


+  (An)2 

360 


C32  34>fJ_1 


264<j)”_2  +  15  9  (})”  _  3 


3  84> 


ft 

N-4 


] 


E3  <|  (An)6  4>VI(S) 


Values  of  [C(f)"3  can  now  be  determined  at  station  N  using 
the  momentum  equation  and  the  extrapolated  values  of  <J>M  ,  (j> '  , 

<J>  • 


C . 2 . 2  Interpolation  Formulas 


(CDn  =  (CDn_x  +  fj[9(C<T)N)E  +  19(00")^.! 


-  5  CC<j.")^_2  +  t  C4>"  )  3  1 


with  an  error  of 

E4  -  "  TTo  (An)5  tcr(0]y 


K  -  «.l  +  -  5*N-2  +  *N-3] 


'N  YN  - 1  24  N 


also  , 
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e5  i  -  tTo  (Ari)5  <t)VI(5> 


and , 


~  ^N-l  +  An<f>’N_i 


(An) 

360 


[38  4>  ” 


+  171^-1  -  36^N-2  + 


7 ^N- 3  ^ 


E6  1  -  OT  (An)6  ^VI(5) 


As  can  be  seen  the  error  in  the  interpolation  formulas 
are  much  less  than  the  error  in  the  extrapolation  formulas. 
Since  the  errors  are  of  opposite  sign,  the  exact  values  must 
be  within  the  bounds  of  the  extrapolated  and  interpolated 
values.  The  solution  can  be  made  more  exact  by  choosing 
smaller  values  of  An. 


o±3 al'o  j  :a  n  '•  a.4  n±  -  -iis  3  ■  >A 


.  .  i  “^or 


77 


O 


U 

z 


FIGURE  1  SCHEMATIC  DIAGRAM  OF  LOW  DENSITY  WIND  TUNNEL 
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FIGURE  2  PHOTOGRAPH  OF  LOW  DENSITY  WIND  TUNNEL 
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FIGURE  3  PHOTOGRAPH  OF  TEST  SECTION  WITH  EQUILIBRIUM  TEMPERATURE  PROBE 
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FIGURE  4  PHOTOGRAPH  OF  TEST  SECTION 
SHOWING  TRAVERSING  MECHANISM 


FIGURE  5  PHOTOGRAPH  OF  IMPACT  PROBE 
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FIGURE  6  SCHEMATIC  DIAGRAM  OF  THE  NOZZLE 


EMF  MILLI  -  VOLTS 
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PRESSURE  (microns) 


FIGURE  7  THERMOCOUPLE  CALIBRATION  CURVE 
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SCALE:  1"  =  3" 


TOP  VIEW 


SCALE  :  1"  =2" 

SIDE  VIEW 


NOTE: 

1.  All  materials  are  aluminum.  The  plates  covering  the  milled 
out  compartments  are  welded  to  the  top  of  the  plate. 

2.  The  inlet  and  outlet  tubes  are  1/4”  aluminium  tubing. 

3.  The  bottom  plate  surface  is  polished  to  a  very  smooth 
finish  and  the  front  of  the  plate  is  machined  to  a  sharp  tip. 


FIGURE  8  SCHEMATIC  DIAGRAM  OF  PLATE  MODEL 
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FIGURE  9  VELOCITY  PROFILE  AT  NOZZLE  EXIT: 


P 


78  MICRONS 


8  5 


FIGURE  10  VELOCITY  PROFILE  AT  NOZZLE  EXIT: 
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FIGURE  11  IMPACT  PROBE  CALIBRATION 


MACH  No.  =  0.64 
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FIGURE  12  SLIP  ALONG  A  COLD  FLAT  PLATE 
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FIGURE  13  SLIP  ALONG  A  FLAT  PLATE: 


T,/T 
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FIGURE  14  TEMPERATURE  JUMP  ALONG  A  FLAT  PLATE: 

=  30.4  MICRONS 
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FIGURE  16  TEMPERATURE  JUMP  ALONG  A  FLAT  PLATE: 
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65.0  MICRONS 


-so 


n/n 


92 


FIGURE  17  VELOCITY  PROFILE  ON  A  COLD  FLA1  PLATE: 

=  0.2  INCHES 
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FIGURE  18  VELOCITY  PROFILE  ON  A  COLD  FLAT  PLATE: 

=  0.4  INCHES 
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FIGURE  19  VELOCITY  PROFILE  ON  A  COLD  FLAT  PLATE: 

=  0.6  INCHES 
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FIGURE  20  VELOCITY  PROFILE  ON  A  COLD  FLAT  PLATE: 

=  0.8  INCHES 
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FIGURE  21  VELOCITY  PROFILE  ON  A  COLD  FLAT  PLATE: 
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FIGURE  22  TEMPERATURE  PROFILE  ON  A  COLD  FLAT  PLATE: 
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FIGURE  23  TEMPERATURE  PROFILE  ON  A  COLD  FLAT  PLATE 

=  0.4  INCHES 
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FIGURE  24  TEMPERATURE  PROFILE  ON  A  COLD  FLAT  PLATE 

-  0.6  INCHES 


x 


100 


8 


i 
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FIGURE  26  TEMPERATURE  PROFILE  ON  A  COLD  FLAT  PLATE: 

=  1.0  INCHES 
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FIGURE  27  VELOCITY  PROFILE  ON  A  FLAT  PLATE: 

P  =  29.0  MICRONS  x  =  0.2  INCHES 
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FIGURE  28  VELOCITY  PROFILE  ON  A  FLAT  PLATE: 

P  =  29.0  MICRONS  x  =  0.4  INCHES 
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FIGURE  29  VELOCITY  PROFILE  ON  A  FLAT  PLATE: 

P  =  29.0  MICRONS  x  =  0.8  INCHES 
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FIGURE  30  VELOCITY  PROFILE  ON  A  FLAT  PLATE: 

=  29.0  MICRONS  x  =  1.0  INCHES 
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FIGURE  31  TEMPERATURE  PROFILE  ON  A  FLAT  PLATE: 

P  =  29.0  MICRONS  x  =  0.2  INCHES 


•  .  .  ■  0 
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FIGURE  32  TEMPERATURE  PROFILE  ON  A  FLAT  PLATE: 

=  29.0  MICRONS  x  =  0.4  INCHES 


P 


108 


FIGURE  33  TEMPERATURE  PROFILE  ON  A  FLAT  PLATE: 

P  =  29.0  MICRONS  x  =  0.8  INCHES 
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FIGURE  34  TEMPERATURE  PROFILE  ON  A  FLAT  PLATE: 

p  =  29.0  MICRONS  x  =  1.0  INCHES 


- 


n/n 


no 


0.8 


0.6 


0.4 


0.2 


0 


/ 

/ 

/ 

/ 

-  / 

/ 

j 

/ 

/ 

f-  o  EXPERIMENT 

. THEORETICAL  WITH  B.C  FROM  EXPERIMENTS 

- - THEORETICAL  WITH  CLASSICAL  B.C. 

I  I  I  I _ I _ I _ I _ I _ I _ L_ . 


MACH  No.  - 

0.64 

P  = 

64.0  microns 

Too 

=  268 °K 

Uco 

=  688  ft./ 

sec. 

X  = 

0.2  in. 

0  0.2  0.4  0.6  0.8 


DISTANCE  (inches) 


FIGURE  35  VELOCITY  PROFILE  ON  A  FLAT  PLATE: 

P  =  6  A . 0  MICRONS  x  =  0.2  INCHES 
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FIGURE  36  VELOCITY  PROFILE  ON  A  FLAT  PLATE: 

=  64.0  MICRONS  x  =  0.4  INCHES 
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FIGURE  37  VELOCITY  PROFILE  ON  A  FLAT  PLATE: 

=  64.0  MICRONS  x  =  0.6  INCHES 
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FIGURE  38  VELOCITY  PROFILE  ON  A  FLAT  PLATE: 

=  6  A . 0  MICRONS  x  =  0.8  INCHES 
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FIGURE  39  VELOCITY  PROFILE  ON  A  FLAT  PLATE: 

=  64.0  MICRONS  x  =  1.0  INCHES 
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FIGURE  40  TEMPERATURE  PROFILE  ON  A  FLAT  PLATE: 

=  64.0  MICRONS  x  =  0.2  INCHES 
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FIGURE  41  TEMPERATURE  PROFILE  ON  A  FLAT  PLATE: 

p  =  64.0  MICRONS  x  =  0.4  INCHES 
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FIGURE  42  TEMPERATURE  PROFILE  ON  A  FLAT  PLATE: 

=  64.0  MICRONS  x  =  0.6  INCHES 
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FIGURE  43  TEMPERATURE  PROFILE  ON  A  FLAT  PLATE: 

=  64.0  MICRONS  x  =  0.8  INCHES 
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FIGURE  44  TEMPERATURE  PROFILE  ON  A  FLAT  PLATE: 

=  64.0  MICRONS  x  =  1.0  INCHES 
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TABLE  1 

EXPERIMENTAL  DATA  AT  THE  WALL  (COOLED  PLATE) 

P  =  64.7  microns  T^  =  267. 9°K 

Uoo  =  689  ft/sec  Tw  =  82.1°K 

T 

JJ_  _ 

Uoo  c  T 

°  K 


Distance  from  Leading  Edge 
inches 


microns 


0.0 

80.5 

0.470 

83.1 

0.1 

77.0 

0.421 

83.1 

0.2 

75 . 2 

0.393 

83.1 

0 . 3 

74.0 

0 .372 

83.1 

0.4 

74.0 

0 .372 

83.1 

0  o  5 

73.7 

0.367 

83.  1 

0.6 

73.  7 

0.367 

83.1 

0.  7 

73.4 

0 . 361 

83.1 

0 . 8 

73.4 

0. 361 

83.1 

0.9 

73.4 

0 . 361 

83.1 

1.0 

73.1 

0.356 

83.1 

- - 
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TABLE  2 


EXPERIMENTAL  DATA  AT  THE 

WALL  (UNCOOLED 

PLATE) 

P  =  30.4  microns 

T 

OO 

=  250.0° 

K 

=  819  ft/sec 

T 

w 

=  301.4° 

K 

Distance  from  Leading  Edge 

inches 

Pt 

microns 

U 

U 

OO 

T 

-1  c 
°K 

u 

T 

X  OO 

0.0 

42.1 

0.903 

302.1 

1.027 

0.1 

37.0 

0.  729 

301.5 

1.078 

0.2 

34.7 

0 .583 

301. 3 

1.114 

0.3 

33.2 

0.475 

301.3 

1.137 

0.4 

32.3 

0 . 392 

301. 3 

1.151 

0.5 

32.3 

0.392 

301.3 

1.151 

0.6 

32.0 

0.279 

301.3 

1.167 

0.7 

31.4 

0 . 322 

301.3 

1.161 

OO 

• 

o 

31.1 

0.279 

301.3 

1.167 

0.9 

31.1 

0.279 

301.3 

1.167 

1.0 

31.1 

0.279 

301.4 

1.167 

. 
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TABLE  3 


EXPERIMENTAL  DATA  AT  THE 

P  =  65  c  0  microns 
Uoo  =  673  ft/sec 

Distance  from  Leading  Edge 
inches 

0.0 
0.1 
0.2 
0. 3 
0.  A 
0.5 
0 . 6 
0.7 
0 . 8 
0.9 
1.0 


WALL  (UNCOOLED  PLATE) 

Too  =  269. 5°K 
Tw  =  300. 5°K 


pt 

microns 

U 

Uoo 

T 

A  c 

°  K 

h. 

T 

A  00 

79.9 

0.893 

301.6 

1.016 

72 . 8 

0.669 

300.8 

1.057 

70 . 1 

0.550 

300.  7 

1.075 

69.5 

0.518 

300 . 6 

1.080 

68.9 

0 .485 

300.6 

1.084 

68.0 

0 . 428 

300.6 

1.091 

66 . 8 

0 . 336 

300 . 6 

1.100 

66 . 8 

0 .336 

300 . 6 

1.100 

66 . 8 

0 . 336 

300 . 6 

1.  100 

66 . 8 

0 . 336 

300 . 6 

1.100 

66.8 

0.336 

300.6 

1.101 
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TABLE  4 


EXPERIMENTAL  DATA  AT  P  =  65.0  MICRONS  (COOLED  PLATE) 


x  =  0.2  inches  T*  =  263  °K 

Uoo  =  664  ft/sec  Tw  =  82.7°K 


Distance  Pt 


inches 


microns 


T-T^ 

Tw“T°° 


.  004 

76.  7 

0.503 

117.  7 

0 . 855 

.  005 

76  .  7 

0.508 

120 . 1 

0.843 

.  006 

76 . 7 

0.509 

120.5 

0. 841 

.  007 

76.7 

0 .511 

121.5 

0 . 836 

.  008 

76.7 

0.515 

123.0 

0 . 828 

.009 

77.0 

0 .522 

124.0 

0 . 824 

.  01 

77.6 

0.536 

125 . 6 

0 . 819 

.  02 

78.2 

0 .576 

139.2 

0 . 751 

.03 

78.5 

0.603 

149.6 

7.000 

.  04 

78.5 

0.629 

162 . 9 

0 ,632 

.05 

78.5 

0.635 

166 . 1 

0.616 

.06 

78.8 

0.653 

172.5 

0 . 585 

.07 

79.0 

0.670 

178.5 

0.556 

.08 

79  3 

0 .683 

182 . 1 

0.539 

.09 

79 . 3 

0.692 

186 . 8 

0 .516 

.  15 

79 . 3 

0 . 706 

194.6 

0.476 

.  20 

79.6 

0.794 

241.  7 

0.240 

.  25 

80 . 2 

0.832 

257.2 

0 . 166 

.  30 

80.5 

0 . 856 

26  7.9 

0.115 

.  40 

81.4 

0.898 

281.6 

0.053 

.  50 

82 . 6 

0.936 

288 . 7 

0 .026 

.  60 

83.2 

0.955 

292 . 7 

0 .011 

.  70 

83.8 

0.971 

295.0 

0.004 

.  80 

84.4 

0.986 

296 . 3 

0 . 002 

.  90 

84.  7 

0.993 

29  7.1 

0.003 

1.00 

85.0 

1 .000 

29  7.6 

0 . 000 

<3  •'  vTAd  JATM3MI  [3^ a? 


;oo. 
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TABLE  5 

EXPERIMENTAL  DATA  AT  P  =  65.0  MICRONS  (COOLED  PLATE) 


x  =  0.4  inches  Too  =  263  °K 

Uoo  =  664  ft/sec  T  =  82.7°K 

w 


Distance  Pt 

inches  microns 


T-Tr 


Tv-1- 


.004 

74.6 

0 . 433 

104.3 

0.917 

.005 

74.6 

0.435 

105 . 1 

0.912 

.  006 

74.6 

0 .437 

105.9 

0 .908 

.  007 

74.6 

0.443 

109.0 

0 . 892 

.  020 

76.1 

0.482 

113.3 

0 . 875 

.030 

76 . 7 

0.535 

113.5 

0 . 773 

.040 

77.0 

0.565 

145.6 

0 .713 

.  050 

77 . 3 

0,574 

147.2 

0 . 706 

.060 

77.6 

0.591 

152.8 

0.678 

.070 

77.6 

0.604 

159 . 5 

0.644 

.080 

77.9 

0 .618 

163.9 

0.623 

.090 

77.9 

0.632 

171.0 

0.586 

.  150 

78  o  8 

0.681 

18  7.5 

0 .507 

.  200 

00 

o 

00 

0 .713 

205 . 8 

0.414 

.250 

79.3 

0 .779 

236  .  7 

0.261 

.  300 

80 . 2 

0 . 822 

250 . 3 

0.199 

.  400 

81.7 

0 . 886 

269 . 2 

0 .114 

.  500 

82  .  3 

0.918 

280.5 

0.062 

.  600 

83.5 

0.955 

287.2 

0 .038 

.  700 

84.1 

0.974 

291.1 

0.023 

.  800 

84.4 

0.984 

293.5 

0.013 

.900 

84.7 

0.992 

295.0 

0.008 

1  o  000 

85.0 

1.000 

295 . 8 

0.006 

• 

£.  >1 

s  .sex 

"  .  x 

* 


oeo. 

oos . 
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TABLE  6 


EXPERIMENTAL  DATA  AT  P  =  65.0  MICRONS  (COOLED  PLATE) 


X 

=  0.6  inches 

T 

00 

263  °  K 

Uoo 

=  664  ft/sec 

T  = 
w 

82 . 7°K 

TT 

T-Too 

Distance 

u 

T  n 

L 

Uoo 

c 

Tw  T°° 

inches 

microns 

°K 

.  004 

74.0 

0 .416 

103.1 

0 . 920 

.005 

74.0 

0.423 

106 . 3 

0 .903 

.006 

74.3 

0.437 

110.1 

0.884 

.007 

74.3 

0.445 

114.1 

0 . 863 

.008 

74.3 

0.446 

114,9 

0.859 

.009 

74,3 

0.452 

117.  7 

0.844 

.010 

74.6 

0.464 

120.  7 

0 . 830 

.020 

75,2 

0.497 

130 . 9 

0.779 

.030 

75.5 

0.517 

13  7.9 

0.  744 

.040 

75 . 8 

0 .534 

143.5 

0 . 716 

.050 

75.8 

0.543 

148.5 

0.690 

*  0  6  0 

76 , 4 

0.565 

153.1 

0.669 

.070 

76.4 

0.576 

159 . 1 

0.639 

.  150 

77.9 

0.633 

172 . 0 

0.580 

.200 

78,5 

0.685 

193.2 

0.475 

.  250 

79 . 3 

0. 776 

234.2 

0.272 

.  300 

80.2 

0 . 812 

243.5 

0.230 

.  400 

81 . 4 

0 . 860 

256 . 2 

0.174 

.500 

82 . 6 

0.910 

2  70.5 

0. 112 

.  600 

83,2 

0.939 

280 . 0 

0 .068 

.  700 

84 . 4 

0.975 

285 . 9 

0.048 

.  800 

84.7 

0.987 

289  .  7 

0 .031 

.900 

85.0 

0.998 

292 . 1 

0.021 

1.000 

85.0 

1.000 

293.5 

0.014 

V  l  i  -.vJ  N 

d00. 

.0 

■ 

■ 

18 
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TABLE  7 

EXPERIMENTAL  DATA  AT  P  =  65.0  MICRONS  (COOLED 


x  =  0.8  inches  T^  =  263  °K 

Uqo  =  664  ft/sec  Tw=  82.7°K 

Dis  tance 
inches 


U 


I Jc 


microns 


c 
°  K 


.004 

00 

o 

CN 

0.394 

102 . 6 

.005 

72  .  8 

0 . 402 

106.5 

.006 

72.8 

0 , 410 

110. 7 

.  007 

72 . 8 

0.415 

113.  7 

.  008 

72,8 

0 . 420 

116 . 6 

.009 

72 . 8 

0,424 

118.  7 

.010 

72.8 

0.427 

120.4 

.  020 

73,1 

0.447 

127.1 

.030 

73.1 

0.451 

129  .  7 

.040 

73,1 

0.463 

136.5 

.050 

74.3 

0.502 

141,9 

.  060 

74.3 

0.509 

145 . 8 

.070 

74 . 3 

0.519 

151.9 

.  100 

7  4.6 

0.561 

172.5 

.  150 

75.5 

0.599 

181.  7 

.  200 

76.7 

0.650 

194.5 

.250 

77.9 

0 . 728 

223,9 

.  300 

79.0 

0.778 

236 . 6 

.  400 

80 . 8 

0 . 836 

246 . 8 

.500 

82.0 

0 .888 

262 . 6 

.  600 

82 . 9 

0 .926 

273.  7 

.  700 

84.4 

0,970 

281.2 

.  800 

84 . 4 

0.979 

286 . 1 

.  900 

84.7 

0.990 

289 . 3 

1.000 

85.0 

1.000 

291.5 

PLATE) 


Tw”T°° 


0 .919 
0 . 899 
0.877 
0.861 
0.846 
0.835 
0 . 826 
0. 792 
0.778 
0.  742 
0 . 720 
0.699 
0.667 
0.562 
0 .519 
0.460 
0 . 316 
0 . 260 
0.220 
0.149 
0 . 100 
0.074 
0.049 
0 .036 
0.027 


a  a= 

. 

‘ 

0^9  .0 

dOO. 

BOO. 

900. 

0X0. 

• 
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TABLE  8 


EXPERIMENTAL  DATA  AT  P  =  65.0  MICRONS  (COOLED  PLATE) 


x  -  1.0  inches 
=  664  ft/sec 


Tm  =  263  °  K 


T  =  82 . 7°K 

w 


Dis  tance 
inches 


microns 


U 


U 


oo 


C 

°K 


T-T 


00 


T  -T 

w  00 


.  004 

71.6 

0. 374 

100 . 8 

0.927 

.  006 

71.6 

0 . 386 

107.8 

0 . 890 

.008 

71.6 

0.394 

112 . 0 

0 . 868 

.010 

71.9 

0.409 

116.4 

0 . 846 

.020 

71.9 

0.427 

126.8 

0. 791 

o  030 

71.9 

0.435 

131 , 5 

0 . 766 

.0  40 

72.2 

0.451 

136.4 

0 . 742 

.050 

72 . 2 

0.457 

140.2 

0  .  722 

.060 

72  .  2 

0.467 

146.6 

0.688 

.  100 

73.4 

0.529 

165.3 

0 .596 

.  150 

74.3 

0.578 

180 . 9 

0 .520 

.  200 

75.5 

0.627 

192.7 

0 . 466 

.250 

76.4 

0.689 

216 . 8 

0.347 

.  300 

77.6 

0 .738 

229  ,  2 

0 . 292 

.  400 

79  .  9 

0.811 

240.3 

0 .252 

.500 

81. 7 

0 .877 

256 . 8 

0 . 181 

.  600 

82 . 3 

0.910 

268.4 

0.127 

.  700 

83.5 

0.949 

276.9 

0.094 

.  800 

83.8 

0 .965 

282.8 

0.067 

.900 

84.4 

0.983 

286  .  7 

0.052 

1.000 

85.0 

1.000 

289.5 

0.042 

. 
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TABLE  9 

EXPERIMENTAL  DATA  AT  P  =  29,0  MICRONS  (UNCOOLED  PLATE) 


x  =0.2  inches 


Too  =  249  °  K 


Uoo  =  872  ft/sec 


Ttt  =  301.  7°K 
w 


Distance  Pfc 

inches  microns 


-T 


00 


.  002 

34.4 

0 .614 

301 . 8 

1. 

,  000 

.003 

34 . 4 

0.614 

301.8 

1. 

,000 

.  004 

34.4 

0 .614 

301.8 

1. 

,  000 

.  005 

34.4 

0.614 

301.8 

1. 

,000 

.  006 

34.4 

0.614 

301.8 

lc 

,  000 

.007 

34.4 

0 .614 

301.8 

1. 

,  000 

.008 

34.8 

0.630 

301.8 

0. 

,  963 

.009 

34 . 8 

0.630 

301.8 

0. 

,963 

.010 

35 . 1 

0.645 

301.8 

0. 

,927 

.  020 

35  .  7 

0.674 

301.8 

0< 

,  857 

.030 

36 . 3 

0 . 702 

301.9 

0, 

,  792 

.040 

37 . 2 

0. 740 

302.0 

0. 

,699 

.050 

37.9 

0 . 764 

302 . 0 

0. 

,635 

.060 

37.9 

0 . 764 

302 . 1 

0. 

,638 

.070 

38.2 

0.775 

302.1 

0. 

,607 

.  080 

38.5 

0 .786 

302.1 

0. 

,5  77 

.090 

38.5 

0 . 786 

302 . 2 

0, 

,579 

.  100 

38.5 

0 .786 

302 . 2 

0, 

,581 

.  150 

39 . 1 

0.808 

302.5 

0  . 

,530 

.  200 

39 . 4 

0 . 818 

302 . 7 

0, 

.  507 

.250 

40.0 

0 .838 

303.0 

0  , 

.  457 

.  300 

40.6 

0.857 

303.1 

0, 

.  405 

.  400 

41.6 

0.885 

303.5 

0, 

.  335 

.  500 

42.8 

0.919 

303. 7 

0  , 

.  241 

.  600 

43.7 

0.942 

303 . 9 

0  , 

.172 

.  700 
.  800 

44.7 

45.3 

0.965 

0.979 

304.1 

304.2 

0, 

0, 

,  106 
.063 

.  900 

45.6 

0.986 

304.2 

0  , 

.043 

1.000 

46.2 

1 .000 

304. 3 

0  , 

.  000 

' 

OM  0 

:  8€ 

v.f  <„  i 

.m 
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TABLE  10 


EXPERIMENTAL  DATA  AT  P  =  29,0  MICRONS  (UNCOOLED  PLATE) 


x  =0.4  inches 


Too  =  249  °  K 


U  oo  =  872  ft/sec 


Tw  =  301 . 7  ° K 


Dis  tance 
inches 


microns 


U 


T-T, 


c 
°  K 


Tw- Tqq 


.  002 

32 . 3 

0 .548 

301 . 6 

1.000 

.003 

32 . 3 

0.548 

301.6 

1.000 

.004 

32.3 

0.548 

301.6 

1.000 

.  005 

32.3 

0.548 

301.6 

1.000 

.006 

32 . 3 

0.548 

301.6 

1.000 

.007 

32 . 3 

0.548 

301.6 

1.000 

.008 

32 . 3 

0.548 

301.6 

1.000 

.009 

32.6 

0.567 

301.6 

0.967 

.010 

32 . 6 

0.567 

301,6 

0.967 

.020 

33.2 

0.601 

301.6 

0 . 903 

.030 

34.1 

0.647 

301 . 6 

0 . 811 

.040 

34.5 

0.661 

301.6 

0.781 

.050 

34.8 

0.675 

301.7 

0.753 

.060 

35.4 

0.701 

301.7 

0.696 

.070 

35  .  7 

0 , 713 

301. 7 

0.668 

.080 

36.0 

0.725 

301.7 

0.641 

.  090 

36 . 3 

0.737 

301. 7 

0.615 

.  100 

36 . 3 

0.737 

301.9 

0 ,618 

.  150 

37.5 

0,782 

302.0 

0.518 

.200 

38.5 

0 . 812 

302 . 2 

0.448 

.250 

39 . 1 

0.831 

302 . 4 

0 . 405 

.  300 

39 . 7 

0 . 850 

302 . 8 

0 . 369 

.  400 

41.3 

0.892 

303.0 

0 . 265 

.500 

42.5 

0.924 

303,4 

0 . 191 

.600 

43.4 

0 .946 

303.6 

0.135 

.  700 

44.0 

0.960 

303.9 

0.103 

.  800 

44.7 

0.974 

304.1 

0 .069 

.900 

45.6 

0.994 

304 . 2 

0 .016 

1.000 

45.9 

1.000 

304.3 

0 . 000 

Sf8  - 


•  e  .0 

w  .tt  , 

> 

voo, 

■ 

Oc  . 

• 

00 A  . 
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TABLE  11 

EXPERIMENTAL  DATA  AT  P  =  29.0  MICRONS  (UNCOOLED  PLATE) 

x  =0.8  inches 
=  872  ft/sec 

Distance  Pfc 

inches  microns 


T^  =  249  °  K 


T  =  301. 7°K 
w 


JJ_ 

u„ 


c 
°  K 


T-T 


OO 


T  -T 

w  00 


.002 
.003 
.004 
.005 
.006 
.007 
.008 
.009 
.010 
.020 
.030 
.040 
.050 
.060 
.070 
.080 
.090 
.  100 
.  150 
.  200 
.250 
.  300 
.  400 
.500 
.600 
.  700 
.  800 
.900 
1.000 


32.6 

0 .482 

32.6 

0.482 

32.6 

0.482 

32.6 

0 .482 

32 . 6 

0.482 

32.6 

0.482 

32.6 

0.482 

32.9 

0.505 

32.9 

0.505 

33.5 

0.546 

33.5 

0.546 

33.8 

0.565 

34.1 

0.584 

34.5 

0.601 

34.8 

0.618 

35.1 

0.634 

35 . 4 

0 .650 

35.4 

0.650 

36 . 3 

0.693 

37.2 

0 .733 

37.9 

0.757 

38.8 

0. 792 

40 . 3 

0.844 

41.9 

0.890 

43.1 

0.924 

44.3 

0.956 

45.0 

0.971 

45.6 

0.986 

46 . 2 

1.000 

301.5 

1.000 

301.5 

1.000 

301.5 

1 . 000 

301.5 

1.000 

301.5 

1.000 

301.5 

1.000 

301.5 

1.000 

301.5 

0.967 

301.5 

0.967 

301.5 

0.903 

301.5 

0.903 

301.6 

0 . 872 

301.6 

0 . 841 

301.6 

0.811 

301.6 

0. 781 

301.6 

0 . 752 

301.6 

0.723 

301.6 

0.723 

301.7 

0.643 

301.  7 

0 .565 

302 . 0 

0.519 

302.1 

0 . 449 

302.5 

0 . 342 

302 . 8 

0.242 

303.1 

0 . 168 

303.5 

0,098 

303.7 

0 . 066 

303.9 

0.033 

304.0 

0 . 000 

"  '  £ 

r  ■  o> 

Q  .  £  0€ 


X  k 

- 


e 

' 

ooe . 
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TABLE  12 

EXPERIMENTAL  DATA  AT  P  =  29.0  MICRONS  (UNCOOLED  PLATE) 


x  =  1.0  inche 


Too  =  249  °  K 


Uoo  =  872  ft/sec 


Tw  =  301 . 7  °  K 


Distance  P 

inches  microns 


T-T 

00 


00 


.002 

32,6 

0.387 

301.4 

1.000 

.003 

32.6 

0.387 

301.4 

1.000 

.004 

32.6 

0.387 

301.4 

1.000 

.005 

32.6 

0.387 

301.4 

1.000 

.006 

32.6 

0.387 

301.4 

1.000 

.007 

32,6 

0.387 

301.4 

1.000 

.008 

32.6 

0.387 

301.4 

1.000 

.009 

32.6 

0.387 

301.4 

1.000 

.010 

33.2 

0.446 

301.4 

0.935 

.020 

33.2 

0.446 

301.4 

0.935 

.030 

33.5 

0.472 

301.4 

0.904 

.040 

33.5 

0.472 

301.4 

0.904 

.050 

33.8 

0.496 

301.4 

0 .873 

.060 

34.1 

0.518 

301.4 

0 . 843 

.070 

34.5 

0.540 

301.4 

0.814 

.080 

34.8 

0.560 

301.5 

0 .785 

.090 

35.1 

0.580 

301.5 

0.756 

.100 

35.1 

0.580 

301.5 

0.756 

.  150 

36.0 

0.633 

301.5 

0.673 

.  200 

36.6 

0.665 

301.6 

0.623 

.  250 

37.5 

0.709 

301.7 

0.549 

.  300 

38.2 

0.737 

301.9 

0 .503 

.  400 

40.0 

0.809 

302 .2 

0 .372 

.  500 

41.3 

0.852 

302.5 

0.293 

.  600 

42.5 

0.892 

302 . 8 

0.219 

.  700 

44.0 

0.936 

303.2 

0 .130 

.  800 

45.0 

0.961 

303.4 

0.080 

.  900 

45.6 

0.977 

303.6 

0.048 

1.000 

46.5 

1.000 

303 . 7 

0 . 000 

' 

■ 

. 
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TABLE  13 


EXPERIMENTAL  DATA  AT  P  =  64.0  MICRONS  (UNCOOLED  PLATE) 


x  =  0.2  inches  T  =  268  °K 

00 

U  =  688  ft/sec  T  =  300. 5°K 

00  w 


Dis  t ance 
inches 


microns 


T-T 


00 


.  002 

70 . 7 

0.608 

300 . 6 

1.000 

.  003 

70 . 4 

0.595 

300 . 6 

1 . 000 

.004 

70.4 

0.595 

300.6 

1.000 

.005 

70 . 7 

0.608 

300.6 

1.000 

.006 

70.7 

0.608 

300 , 6 

1.000 

.007 

70 . 7 

0.608 

300.6 

1 . 000 

.008 

70 . 7 

0 . 608 

300.6 

1.000 

.  009 

71.0 

0 .620 

300 . 6 

0.970 

.010 

71.3 

0.632 

300.6 

0.940 

.020 

72.2 

0.667 

300.6 

0 . 853 

.030 

72.8 

0.689 

300.6 

0.797 

.040 

73.1 

0 . 700 

300 . 8 

0 .778 

.050 

73.4 

0 .710 

300 . 8 

0 . 752 

.060 

74.0 

0 ,731 

300.9 

0.699 

.070 

74.3 

0 . 740 

300.9 

0.672 

.080 

74.6 

0 . 750 

301.0 

0  o  650 

.090 

74.6 

0 .750 

301 . 1 

0.655 

.  100 

74.9 

0 . 760 

301.2 

0 .633 

.  150 

76.1 

0.796 

301.6 

0.548 

.  200 

76  .  7 

0.814 

301.9 

0.510 

.  250 

77  .  3 

0 . 831 

302 . 2 

0.473 

.  300 

77.6 

0 . 840 

302 . 8 

0.478 

.  400 

79.3 

0 . 886 

303 . 3 

0.353 

.  500 

80.2 

0.908 

303.5 

0.293 

.  600 

81.4 

0.936 

303 . 7 

0.207 

.  700 

82 .0 

0.950 

303.9 

0 .171 

.  800 

83.5 

0.982 

303.9 

0.063 

.900 

83.8 

0.990 

304.0 

0.043 

1 . 000 

84.4 

1.000 

304.0 

0.000 

0::  J)  3 WO  M  0.  d  •  S'  TA  ATACI  .  ‘  TWSTM 1  Si  'J?: 

' 


80*. 0 

V  .0^ 

t  .0^ 

. 

010. 

s.  sv 

■ 

V 

d.  AV 

• 

OOd. 

oos;. 

0^9.0 
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TABLE  14 


EXPERIMENTAL 

DATA  AT 

P  =  64.0 

MICRONS 

(UNCOOLED  PLATE) 

X 

=  0.4 

inches 

Too  = 

268  °  K 

Uoc 

,  =  688 

ft/ sec 

T  = 
xw 

300 . 5°K 

Distance 

inches 

P 

U 

T 

Y  t 

microns 

Uoo 

1  C 

°  K 

Tw-Too 

.002 

69 . 8 

0.559 

300 . 5 

1.000 

.003 

69.5 

0.545 

300.5 

1.028 

.  004 

69.5 

0.545 

300.5 

1.028 

.005 

69 . 8 

0.559 

300 . 5 

1.000 

,006 

69 . 8 

0.559 

300.5 

1.000 

.007 

69  .  8 

0.559 

300.5 

1.000 

.008 

69.8 

0.559 

300.5 

1 . 000 

.009 

69.8 

0.559 

300.5 

1 . 000 

.010 

70 . 1 

0.573 

300.5 

0.972 

.  020 

70.4 

0.586 

300.5 

0 .947 

.030 

71.0 

0.612 

300.6 

0.893 

.040 

71.3 

0.625 

300 . 6 

0.866 

.050 

71.9 

0.648 

300.6 

0 .816 

.  060 

72.5 

0.671 

300  .  7 

0 . 766 

.070 

72.5 

0.671 

300 . 7 

0 . 766 

.080 

73.1 

0,693 

300 . 7 

0 .718 

.090 

73.1 

0.693 

300 . 8 

0 . 722 

.  100 

73.4 

0 . 703 

300 . 8 

0.697 

.  150 

74.9 

0.753 

301.2 

0.589 

.200 

76 . 1 

0 . 789 

301.6 

0.510 

.250 

77 . 3 

0 . 824 

301.9 

0 .430 

.  300 

78.5 

0.856 

302 .1 

0 . 350 

.  400 

79.6 

0 .887 

303.0 

0 . 300 

.500 

80 . 8 

0.916 

303.5 

0.237 

.600 

82 .0 

0.943 

303 . 8 

0 . 166 

.  700 

82.9 

0.963 

303.9 

0 .113 

.  800 

83.5 

0.976 

304.0 

0.074 

.900 

84.1 

0.988 

304.1 

0.039 

1.000 

84.7 

1.000 

304 . 1 

0 .000 

.  x  . 

GOO. 

■  .'L. 

e .  ooe 

. 

• 

)  COG 

oeo. 

. 

OOG. 

o.i<n 

ooe. 
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TABLE  15 


EXPERIMENTAL 

DATA  AT 

P  =  64.0 

MICRONS 

(UNCOOLED  PLATE) 

X 

=  0.6 

inches 

Too  = 

268  °K 

Uoc 

,  =  688 

f  t  /  s  ec 

T  = 
xw 

300. 5°K 

Dis  tance 

p 

U 

T 

T-T 

X  x  00 

Uoo 

T  C 

8 

H 

1 

H 

inches 

microns 

°  K 

.  002 

6  8.6 

0 . 495 

300 . 6 

1.000 

.003 

68.6 

0.495 

300 . 6 

1.000 

.004 

6  8.6 

0.495 

300 . 6 

1.000 

.005 

68.6 

0.495 

300 , 6 

1.000 

.006 

6  8.6 

0 . 495 

300.6 

1.000 

.007 

68.3 

0.478 

300.6 

1.027 

.008 

68.9 

0.512 

300 . 6 

0.974 

.009 

68.9 

0.512 

300.6 

0.974 

.010 

69.2 

0.527 

300.6 

0.948 

.020 

69.5 

0,542 

300.6 

0.922 

.030 

69.8 

0.557 

300.6 

0 . 896 

.040 

70.7 

0.598 

300.6 

0 . 820 

.050 

71.0 

0.611 

300.6 

0 . 798 

.060 

71 . 3 

0.624 

300.6 

0.774 

.070 

71.3 

0.624 

300 . 7 

0.775 

.080 

71.3 

0.624 

300 . 7 

0,775 

.090 

71.6 

0.636 

300 . 7 

0 . 754 

.100 

71.6 

0.636 

300 . 8 

0 .755 

.  150 

74 . 3 

0 ,735 

301.0 

0.555 

.200 

75.2 

0 . 764 

301 . 3 

0.501 

.250 

76.4 

0.801 

301 . 7 

0,428 

.  300 

77.9 

0 . 843 

301.9 

0.331 

.  400 

79.3 

0.883 

302.4 

0.250 

.500 

80 . 8 

0.920 

303.2 

0.186 

.600 

82 .0 

0.948 

303.6 

0.126 

.  700 

82 .6 

0.962 

303.9 

0.099 

.  800 

83.8 

0.988 

304.1 

0.033 

.900 

84.1 

0.994 

304.1 

0.018 

1.000 

84.4 

1,000 

304 . 1 

0 , 000 

e  :  Hi a AT 


' 


.  a-  N 

. 

■ 
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TABLE  16 


EXPERIMENTAL  DATA  AT  P  =  64.0  MICRONS  (UNCOOLED  PLATE) 


x  =0.8  inches 
U oo  =  688  ft/ sec 


Too  =  26  8  °K 


Tw  =  300. 5°K 


Distance 

inches 


microns 


U 


x  c 
°K 


T-T 


oo 


T  -T 
iw  00 


.  002 

67.4 

0 . 422 

300.6 

1.000 

.  003 

67.1 

0.401 

300 . 6 

1.025 

.004 

67.1 

0 . 401 

300.6 

1.025 

.005 

67.1 

0 . 401 

300.6 

1.025 

.006 

67.1 

0.401 

300.6 

1.025 

.007 

67.1 

0.401 

300.6 

1.025 

.008 

67.1 

0.401 

300.6 

1.025 

.009 

67.7 

0.442 

300.6 

0.975 

.010 

67.4 

0.422 

300.6 

1.000 

.020 

6  8.3 

0.478 

300.6 

0.928 

.030 

68.6 

0.495 

300 . 6 

0.904 

.040 

68.9 

0.512 

300.6 

0 . 880 

.050 

68.9 

0.512 

300.6 

0 .881 

o  060 

69  .  2 

0.527 

300 . 6 

0.857 

.070 

69.5 

0.542 

300.6 

0 .834 

.080 

69.8 

0.557 

300 . 7 

0.811 

.090 

70 . 4 

0.585 

300 . 7 

0 . 766 

.  100 

71 . 3 

0.624 

300 . 8 

0  .  702 

.150 

71.9 

0.649 

301.0 

0.666 

.200 

73.1 

0.694 

301 . 1 

0.588 

.250 

74.3 

0 .736 

301.5 

0.518 

.300 

76.4 

0 . 801 

301.8 

0 . 389 

.  400 

77.9 

0.844 

302.2 

0 .312 

.  500 

79.6 

0.891 

302 .9 

0.228 

.600 

81.1 

0.928 

303.4 

0.158 

.  700 

82  .  3 

0.955 

303.8 

0 . 102 

.  800 

83.5 

0.981 

304 .0 

0.042 

.900 

84 . 1 

0.994 

304.1 

0.016 

1.000 

84.4 

1.000 

304 . 1 

0 .000 

a.ooe 

- 4  a 
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d.  00 i 
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TABLE  17 

EXPERIMENTAL  DATA  AT  P  =  64.0  MICRONS  (UNCOOLED  PLATE) 


x  =  1.0  inches  T^  =  268  °K 

Uoo  =  688  ft/sec  Tw  =  300. 5°K 


Distance  P, 

t 


inches 


microns 


U 

Uoo 


T 

c 


°K 


T-T 


OO 


-T 


OO 


,002 

65 . 3 

0.270 

300 . 6 

1.000 

.  003 

65 . 3 

0.270 

300 . 6 

1 . 000 

.  004 

65.3 

0.270 

300.6 

1.000 

.005 

65 . 3 

0.270 

300.6 

1.000 

.006 

65  .  3 

0.270 

300.6 

1.000 

.007 

65 . 3 

0.270 

300 . 6 

1.000 

.  008 

65.9 

0 . 328 

300.6 

0 . 956 

.009 

65 . 9 

0.328 

300 . 6 

0.956 

.010 

65 . 3 

0.270 

300.6 

1.000 

.  020 

65.9 

0.328 

300 . 6 

0.956 

.030 

6  6.2 

0.353 

300 . 6 

0.934 

.040 

67.4 

0.439 

300 . 6 

0.849 

.050 

67.4 

0.439 

300 , 6 

0.849 

.060 

67.7 

0 . 458 

300 . 6 

0.828 

.070 

68.0 

0 .475 

300 . 6 

0 . 808 

.080 

68.6 

0 .509 

300 . 6 

0.767 

.090 

6  8.6 

0.509 

300 . 7 

0  768 

.  100 

6  8.6 

0.509 

300 . 7 

0 . 768 

.  150 

70 . 7 

0 . 608 

300 . 9 

0 .638 

.200 

71,6 

0 .645 

301 . 1 

0 .588 

.250 

72.5 

0.679 

301 . 3 

0.542 

.  300 

73.4 

0 .711 

301 . 5 

0 . 495 

.  400 

75 . 8 

0  .  788 

302 . 0 

0 .373 

.  500 

76  .  7 

0 . 815 

302 . 6 

0.343 

.600 

79 . 3 

0 . 886 

303.1 

0.219 

.  700 

81.4 

0 .936 

303.6 

0 . 128 

.  800 

83.8 

0 .988 

303.8 

0.021 

.900 

84.1 

0 .994 

304.0 

0.014 

1.000 

84.4 

1.000 

304.0 

0.000 
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